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Zusammenfassung 

Gegenstand der Arbeit ist die numerische Untersuchung einer iiber das Spek- 
tralwirkungsprinzip definierten nichtkommutativen Feldtheorie. 
Ausgangspunkt dieser Konstruktion ist ein (als harmonisch bezeichnetes) spek- 
trales Tripel (A4, H4, D4). Dabei ist A4 die 4-dimensionale nichtkommuta- 
tive Moyal- Algebra und D 4 ein selbstadjungierter (Dirac-) Operator auf dem 
Hilbert-Raum "H 4 , so dass V\ der Schrodinger-Operator des 4-dimensionalen 
harmonischen Oszillators wird. Die Konstruktion dieser Daten basiert auf 
einer 8-dimensionalen Clifford-Algebra. Fiir das Produkt aus dem Tripel 
(A4, "Hzi, D4) mit einem matrixwertigen spektralen Tripel wird analog zum 
Standardverfahren der nichtkommutativen Geometrie die Spektralwirkung 
berechnet. 

Die Renormierungstheorie assoziiert zur Spektralwirkung ein (z.T. nur for- 
males) Wahrscheinlichkeitsmafi, deren zugehorige Korrelationsfunktionen eine 
Feldtheorie definieren. In der storungstheoretischen Variante wird das 
Wahrscheinlichkeitsmafi als formale Potenzreihe konstruiert. Voraussetzung 
dafiir ist die explizite Kenntnis der Losungen der Euler-Lagrange-Gleichungen 
zur Spektralwirkung. Fiir das betrachtete Modell erweist es sich als unmoglich, 
diese Losungen (als Vakuum bezeichnet) zu gewinnen. 

Ein alternatives Verfahren besteht in der Diskretisierung aller Variablen 
und der numerischen Untersuchung des Verhaltens der Korrelationsfunktio- 
nen bei Verfeinerung der Diskretisierung. Die Diskretisierung des Modells 
wird durch die Matrix-Basis des Moyal-Raums und Beschrankung auf endliche 
Matrizen erreicht. Durch Monte-Carlo-Simulationen werden wichtige Korre- 
lationsfunktionen wie die Energiedichte, die spezifische Warme sowie einige 
Ordnungsparameter untersucht, sowohl in Abhangigkeit von der Grofie der 
Matrizen als auch von den unabhangigen Parametern des Modells. Dabei 
werden trotz der grofien Komplexitat der approximierten Spektralwirkung 
verlafiliche numerische Resultate erzielt, die zeigen, dafi eine numerische Be- 
handlung dieser Art von Modellen in der Matrix- Moyal-Basis moglich ist. 
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Abstract 

The object of this work is the numerical investigation of a non-commutative 
field theory defined via the spectral action principle. Starting point of this 
construction is a spectral triple (A4, H4, T>±) referred to as harmonic. Here, 
v4.4 is the 4-dimensional (noncommutative) Moyal algebra, and V4 is a selfad- 
joint (Dirac-)operator on the Hilbert space H4 such that V\ is the Schrodinger 
operator of the 4-dimensional harmonic oscillator. The construction of these 
data relies on an 8-dimensional Clifford algebra. In analogy to the standard 
procedure of non-commutative geometry, the spectral action is computed for 
the product of the triple (Aa, H4, T>^) with a matrix-valued spectral triple. 

Renormalization theory associates to the spectral action a (often only for- 
mal) probability measure. Its associated correlation functions define then a 
field theory. In the perturbative approach this measure is constructed as a 
formal power series. This requires explicit knowledge of the solutions of the 
Euler-Lagrange equations for the spectral action. For the model under con- 
sideration, it turns out impossible to obtain these solutions. 

An alternative approach consists in a discretization of all variables and 
a numerical investigation of the behavior of the correlation functions when 
the discretization becomes finer. For the model under consideration, the dis- 
cretization is achieved in the matrix basis of the Moyal algebra restricted to 
finite matrices. By Monte Carlo simulation we study several important corre- 
lation functions such as the energy density, the specific heat and some order 
parameters, in dependence both of the matrix size and of the independent pa- 
rameters of the model. Despite the complexity of the approximated spectral 
action, some reliable numerical results are obtained, showing that a numerical 
treatment of this kind of models in the Moyal matrix basis is possible. 
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Introduction 



The main object of this work is a particular non-commutative field theory 
which is derived using the spectral action principle and then treated numeri- 
cally. Non-commutativity can be found in many fields of physics like quantum 
field theories, string theory [5], condensed matter physics. The first application 
of non-commutativity into physicss is dated from the middle of the last century 
inspired by the ideas of quantum mechanics, where starting from classical me- 
chanics, the commutative algebra of functions on the phase space is replaced 
by a non-commutative operator algebra on a Hilbert space. The duality be- 
tween ordinary spaces M and proper commutative algebras is expressed by the 
Gel'fand-Naimark theorem which states the fact that the algebra of all con- 
tinuous functions on M is the only possible type of commutative C*-algebra. 
Additionally, given a commutative C*-algebra C, it is possible to reconstruct a 
Hausdorff topological space M in order to obtain that C is the algebra of con- 
tinuous functions on M. The study of commutative C*-algebras is equivalent 
to the study of topological Hausdorff spaces. The previous duality has inspired 
the identification, in non- commutative geometry, of some algebraical objects 
as a category of non-commutative topological spaces. Alain Connes |31J, one 
of the founders of non-commutative geometry, has proposed a candidate for 
the objects of such category, the spectral triples (A, %, V) [32], composed by 



an algebra A, an Hilbert space H on which A is represented and an self adjoint 
operator T>. In fact, every compact oriented Riemannian manifold can be used 
to define a spectral triple, this kind of manifold M characterizes a symmetric 
Dirac type operators on self-adjoint Clifford module bundles over M. Connes, 
after a conjecture in 1996 [36] and some considerable attempts of Rennie and 
Varilly [67], proved the so called reconstruction theorem [33] for commutative 
spectral triples satisfying various axioms, showing that exists a compact ori- 
ented smooth manifold X such that A = C°°(X) is the algebra of smooth 
functions on X and every compact oriented smooth manifold emerges in this 
way. Pushed by the aim of reformulating the standard model of particles in a 
non-commutative way [351 [36], Connes has introduced the almost-commutative 
spectral triple extending the axioms of the reconstruction theorem to a non- 
commutative algebra. The almost-commutative spectral triples are defined 
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as the non-commutative Cartesian product of a commutative spectral triple 
of a compact spin manifold, with a spectral triple where the Hilbert space is 
finite-dimensional, this triples are often labeled as finite spectral triple. 

A field theory can be interpreted [3] as a theory concerning maps <fi (usu- 
ally referred as fields) between the space-time M' and a target space M. On 
this spaces are defined some structures, they can be Riemannian or pseudo- 
Riemannian manifold depending on the signature of the space-time. A funda- 
mental role in the field theory is played by the action S(4>) which is a functional 
of the elementary fields in the theory. In the classical theory the aim is to study 
the extrema of the action functional 5S((j)o) = in order to obtain the solu- 
tion of the equations of motion. Quoting the intuitive Feynman formulation of 
quantum field theory [61] , the action is used to study the functional integrals: 



where F((f>) is a functional of </> and Dtp is the functional measure. The previous 
expression in general is not well defined and should be considered as an ap- 
proximate expression. However, the functional integral approach is very useful 
in studying the quantum field theory connection with the expectation values 
in statistical mechanics; all field configurations contributes in the estimation 
of the previous functional each whit probability amplitude P((f>) ~ e~n s ^\ 
More formally the action functional is used to compute the correlation func- 
tions which are one the main physical output of a quantum field theory with 
a close connection to the experimental measurements. In the path integral 
formalism [62] the correlation function of the fields <pi(xi) ■ ■ ■ 4> n {x n ) is given, 
for the euclidean case, by the functional integral: 



where is the path integral measure _D$ on the space of configuration of the 
elementary fields. In the Lorentz signature the last term is replaced by e lS . A 
crucial point in the quantum field theory is the formalization of the previous 
integral in order to obtain all correlation functions of the theory well defined. 
This procedure is called renormalization and is achieved using a perturbative 
approach. A safe perturbative analysis can only be conducted after expanding 
the action around its vacuum, of course this requirement needs the explicit 
expression for the vacuum usually obtained minimizing the action. It is clear 
that for any quantum field theory the explicit determination of vacuum is a 
indispensable step before the perturbative study of its renormalizability can 
be done. 

In order to define an action for a non- commutative geometry theory Connes 
and Chamsedinne introduced a general formalism for spectral triples, the spec- 
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tral action principle |53j. The term spectral come from the fact that it depends 
only on the spectrum of the Dirac operator and it takes the form 

Tr( X (V/A)) 

where A is a real parameter and it fixes the energy scale. The x is a differen- 
tiable function x '■ ^+ — ► ^+ of sufficiently fast decrease such that the spectral 
action converges. 

Using this approach the standard model of particles has been reformulated, 
including a Riemannian formulation of gravity. The algebra of the spectral 
triple used is defined by the tensor product of C°°(M), the regular functions 
on a manifold M, times a matrix algebra of finite dimension. Formally the 
natural group of invariance of the standard model, including gravity, is the 
semidirect product: 

G = U x. Diff(M) 

Were U = C°°(M, U{1) x SU(2) x 377(3)) is the group of local gauge trans- 
formations. The total group of invariance admits U as a normal subgroup. 
A.Connes has showed, defining A = C°°(M) <g> (C © H © M 3 (C)), that the 
automorphisms group Aut(.A) of the non-commutative algebra A admits the 
inner group Int(^4) as a normal subgroup and Aut(A) ~ G. The algebra 
Af = (C ©H© M 3 (C)) is finite and H is the algebra of quaternions. The rest 
of the spectral geometry is defined by the action of A on % and by suitable 
Dirac operator operator T>. In [53J is showed that the total spectral triple is 
given by: 

(A, H, V) = (C°°(M), L 2 (M, S), $ M ) © (A F , H F , V F ) 

with 

U = L 2 (M, S) ®H F , V = $ M © 1 + 75 © V F 

The algebra A F is finite dimensional so the dimension of the corresponding 
Hilbert space H F must be finite dimensional [51]. The final step, in order 
to obtain a field theory, is to define an action, this task is achieved using 
the spectral action principle and the invariance under the symmetry group 
is implement fluctuating the Dirac operator. The subgroup Int(^4) of inner 
automorphisms is a normal subgroup and the group Aut (A) of diffeomorphisms 
falls in equivalence classes under Int(y4.). This induces a natural foliation into 
equivalence classes in the space of metrics. The internal fluctuations [38J of a 
given metric are given by the formula: 

V = V + A + JAJ~\ A = J2 a A D ^ b *} a i> he A 

where J is the real structure. In this way starting from (A,H,'D ), instead 
modify the representation of A in H, it is changed the operator V where A 
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is a self-adjoint operator in % of the form A = ^cijfDo,^]- The spectral 
action principle applied to inner fluctuations reproduces the bosonic part of 
the model, the gravity is naturally present in the model, while the other in- 
teractions are encoded in the matrix algebra of the total spectral triple. The 
non-commutativity of matrices corresponds to the non-abelianity of the gauge 
theory defining a so called Yang-Mills theory. 

The previous example of application of non-commutative geometry to a 
field theory is founded on the use of a spectral triple where the algebra is 
almost-commutative, a very interesting task is to formulate a field theory for 
a truly non-commutative algebra. A first attempt was obtained replacing in 
the usual field theory action the point-wise multiplication of the fields with 
a non- commutative one, namely a ^-product. The fields now belongs to Rq, 
a vector space defined by an enough regular class functions on IRq equipped 
with the Moyal product: 



Where is a skew-symmetric matrix. Unfortunately all the attempts to renor- 
malize such quantum field theories on the non-commutative M 4 failed and 
these models show a phenomenon called UV/IR- mixing [T3l [TH [15]. A great 
step towards the non-commutative field theory was made when H.Grosse and 
R.Wulkenhaar jS], found a non-commutative y9 4 -theory renormalizable action 
which develops additional marginal coupling, corresponding to an harmonic 
oscillator potential for the real-valued free field <p on Wq : 



Where x = 20 _1 • x, A G M. Q G [0, 1], and fi is a real parameter. Using 
the Moyal matrix base, which turns the ^-product into a standard (infinite) 
matrix product, H.Grosse and R.Wulkenhaar were able to prove the perturba- 
tive renormalizability of the theory [13]. Afterward, R.Wulkenhaar et al. [18] 
found an alternative simpler normalization proof using multi-scale analysis in 
matrix base, showing the equivalence of various renormalization schemes. A 
last, but useful, renormalization proof was formulated using Symanzik type 
hyperbolic polynomials J2TJ. It is worthwhile to mention [2T| [22]. 

The non- commutative model treated in this thesis is a sort of extension, via 
spectral action principle, of the scalar W-G model, in which we are interested 
to formulate a Yang-Mills theory in renormalizable way on Moyal space. From 
the previous discussions we can expect that usual Yang-Mills theory on Moyal 
space without modifications of the action by something similar to an oscillator 
potential, to be not renormalizable [16]. Additionally, the Moyal space with 
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usual Dirac operator is a spectral triple, the corresponding spectral action was 
computed in [53], with the result that it is the usual not renormalizable ac- 
tion on Moyal plane. In jUJ H.Grosse and R.Wulkenhaar, in order to obtain 
a gauge theory with an oscillator potential via the spectral action principle, 
used a Dirac operator constructed using the statement T> 2 = H where the four 
dimensional Laplacian is substituted by the four dimensional oscillator Hamil- 
tonian H = —A + f2 2 ||a;|| 2 . The idea behind is that the spectral dimension is 
defined through the Dirac operator so the spectral dimension defined by such 
Dirac operator is related to the harmonic oscillator phase space dimension. 
It turns out that to write down an Dirac operator, so that its square equals 
the 4D harmonic oscillator Hamiltonian, is an easy task using eight dimension 
Clifford algebra. In addition, can be shown that using this Dirac operator on 
4D- Moyal space, is possible define an eight-dimensional spectral triple. In this 
thesis is used the approach described in [57] in which the Dirac operator is 
constructed using n-dimensional bosonic and fermionic creation and annihila- 
tion operators. After defined the Dirac operator with the desired spectrum it 
is considered the total spectral triple as the tensor product of the "oscillat- 
ing" spectral triple (Aj, H^, T>^) with an almost-commutative triple and then 
is perform the previous described procedure of non- commutative geometry 
to compute the spectral action. We notice that matrix algebra introduces 
an extension of the standard potential in the commutative case, in fact the 
scalar field (j> an d the fields are present together in a potential of the forrrfl 
(aX^ * X^ + /3<p * ip - l) 2 , with a, f3 G R + and X^x) = (O^ 1 )^ + A^x) is 
a covariant coordinate. These two additional terms, the integral over X^ * X M 
and over its square, were conjectured in [16]. 

The high non-triviality of the vacuum makes very difficult to explicit the 
vacuum configuration of the system in [7_T] A. de Goursac, J.C Wallet, and 
R. Wulkenhaar, using the matrix base formalism, have found an expressions 
from vacuum solutions deriving them from the relevant solutions of the equa- 
tions of motion. Although, the complexity of the vacuum configuration makes 
the perturbative approach very complicated, in order to conduct some inves- 
tigations in this thesis will be considered a non-perturbative scheme using a 
discretized matrix model of the action in which the fields become matrices, 
the star product become the matrix multiplication and the integral turns in a 
matrix trace. 

Now comes in to play the numerical treatment, the standard method is to 
approximate the space by discrete points, for example using a lattice approxi- 
mation and then calculate the observables over that set of points [76]. Since an 
approximation in the position space is not suitable due to the oscillator factor 
of the Moyal product, instead the lattice approximation, will be used the ma- 

1 Einstein notation on repeated indices is used. 
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trix Moyal base, which was already used in the first renormalization proof of 
<^ 4 -model restricted to finite matrices. Hence, will be performed a Monte Carlo 
simulation studying some statistical quantity such the energy density and spe- 
cific heat varying the parameters Q, and gathering some informations 
on the various contributions of the fields to the action. The simulations are 
quite cumbersome due the complexity of the action and the number of inde- 
pendent matrices to handle but using particular algorithm we are able to get 
an acceptable balance between the computation precision and the computa- 
tion time. For the simulations is applied a standard Metropolis-Monte Carlo 
algorithm [75] with various estimators for the error (see appendix B) and for 
the autocorrelation time of the samples. The initial conditions of the Markov 
chain are chosen randomly in the phase space, they can be of two types; hot 
initial conditions, which are configurations far from the minimum, and cold 
start conditions, which correspond to configurations close to the minimum of 
the action. In general we chose the range of parameters in order to avoid prob- 
lems with the thermalization process, obtaining numerical simulations where 
is enough to wait a relative small number of Monte Carlo steps to compute 
independent results from the initial conditions. For the model under study, we 
are interested on the continuous limit that correspond to matrices of infinite 
size. We will consider various size of the matrices expecting a stabilization of 
the values of observables like the energy density, increasing the matrix size. 
In order to find same possible phase transitions will be used the specific heat 
which is a measure of the dispersion of the energy. The phase transitions are 
registered as peaks of the specific heat, increasing the matrices size. Beside, 
we define others quantities, such as (Tr[yj] 2 ) and (5^ n =o Iv^wl 2 ) which are used 
as order parameters [81], we have also defined their susceptibilities but are 
not collected due the high number of samples required to obtain a sufficient 
precision 



Plan of the work 

This work is structured in order to introduce the reader to the model treated 
and to the numerical analysis, using an heuristic approach. The first chapter 
and the third one are introductory chapters in which are explained respectively 
the fundamentals non-commutative geometry and the main concepts of Monte 
Carlo numerical simulation. In particular in the first chapter we will introduce 
the mathematical tools and the machinery required to formulate our spectral 
model. In the first part of chapter one will be described the basics of non- 
commutative geometry, Weyl-Wigner map. Further will be introduced some 
notions of spectral geometry and in the end will presented a simple example of 
application of the spectral action principle. The second chapter is devoted to 
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construct the W-G spectral action model. Will be described how to construct 
the harmonic Dirac operator starting from the one-dimensional case and then 
generalizing it using annihilations and creations operators for the bosonic and 
fermionic sectors. Having the 4D-Dirac harmonic operator with harmonic 
oscillator spectrum and extending it to an spectral triple, will be considered 
the tensor product of the non- commutative triple with a finite Connes-Lott 
type spectral triple [37J. Following the previous described standard procedure 
of non-commutative geometry, to obtain a "gauged" Dirac operator, we will 
fluctuate the total Dirac operator, thus we will proceed to compute the spectral 
action in which are present two U(l)-Moyal Yang-Mills fields unified with a 
complex Higgs field. In the end of the chapter will be discussed some aspects 
of the resulting action in particular of vacuum and the needs for numerical 
treatment. The third chapter is an introduction to Monte Carlo analysis, will 
be set up all the ingredients required to conduct a numerical simulation of 
the model treated in the previous chapter. Will be explained the basics of 
the Monte Carlo simulation, focusing ourself on the application in the field 
theory. In the first sections of chapter three will be briefly discussed the path 
integral formulation, this formulation is essential to connect the field theories 
to statistical systems in which the Monte Carlo methods are born. Then will 
be introduced the Metropolis algorithm used to produce a Markov chain. In 
order to resume all the previous concepts and to show an example of phase 
analysis will be presented an application of numerical simulation on the well 
know Ising model. Right before the presentation of the numerical results, in 
the chapter 4, will be showed how the Monte Carlo it is implemented in our 
case starting from the discretization scheme introducing the Moyal base. In the 
second part, in order to define the observables of the upcoming Monte Carlo 
simulation, will be designated the expectation values, statistical quantities like 
energy, specific heat and the order parameters. Finally, in the last chapter will 
be showed the numerical results for the previous defined quantities for the full 
4-dimensional model and for the 2-dimensional one. 
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Chapter 1 

Introduction to Connes-Lott 
models 

In this chapter we introduce some mathematical tools taken from non-com- 
mutative geometry in order to define a spectral triple, a spectral action and 
Connes-Lott models. In the first part of this chapter we describe the fun- 
damentals of non- commutative geometry in particular the Weyl-Wigner map, 
further will be introduced some notions of spectral geometry. 

1.1 Non-commutative geometry and 
Weyl-Wigner map 

In this section we will discuss, as an example of non-commutative space, the 
fate of the classical phase space under the process of quantization. This is 
a very large topic and we limit ourselves to some aspects. In particular we 
wish to introduce a way to quantize the phase space (and in general an L 2 (M. n ) 
space) in which the quantized space can be seen as either a set of operators on 
an Hilbert space or as a deformation of the product of functions on the classical 
space. At least in the simple cases, the non-commutative algebra of a quantum 
phase space can be taken as the one generated by the position and momentum 
operators acting on a separable Hilbert space. We now establish a connection 
between the classical and quantum phase spaces, at least for systems with 
a well defined classical counterpart. We want to associate to each classical 
observable an bounded operator on H. This will be done by the Weyl map pQ. 
Starting from the correspondence principle, an immediate problem is to solve 
an ordering ambiguity, beside p and q are unbounded operators, additionally 
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arises problems about the definition of their domains. To solve these problems 
Weyl has suggested to introduce a base in the operator space on L 2 (IR) using 
the unitary operators: 

U(q)=e* qP (1.1) 
V(p) = e^ pQ (1.2) 

Where P and Q are respectively the momentum operator and the position one. 
Being bounded, their domains can be extended to the whole Hilbert space 
L 2 (R). In general it is convenient to use the base defined by the following 
unitary operators: 

W(a,P) = e ^ aP+ ^ (1.3) 

Where a and are two real parameters. A generic operator using this base 
can be expanded like: 

F(P,Q) = -^-[ g(a,(3)e^ p+ ^dad(3 (1.4) 
2-Kh J M 2 

Beside this expansion let us considering the base of L 2 (R 2 ): 

w{a,P) = e^ ap+ ? q) (1.5) 
using this base any function can be written as 

f(p,q) = 7^r[ f(a,(3)e^ ap+ ^dadp (1.6) 
2nh J R2 

where f(a, (3) is the Fourier transform of f(p, q). Now, thanks to the expansion 
(jl.6p and (II. 4p , it is possible to build an application which maps each function 
in L 2 (R 2 ) to an bounded operator. In order to construct this map we use the 
function f(a, ft) in the operatorial expansion (II. 4p . In such a way the ordering 
problem is solved, in fact the association is no more ambiguous. 

Formally the Weyl map is the application: 

Q : L 2 (R 2 ) ->• B{U) 

defined by 

f(p,q)^S}(f)(P,Q) = -L [ f(a,(3)W(a,f3)dadf3 (1.7) 

2?m J R 2 

with 

iW) = A/ f(p,q)e-^ ap+ P q) dpdq (1.8) 
2-Kh J R 2 
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This particular Weyl map selects a certain ordering: the symmetric order- 
ing or the so called Weyl's ordering. Other ordering are possible, for example 
normal ordering or Wick ordering, introducing into (ll.7p a weight function [6J . 
In addition, Weyl map is linear and satisfies the condition: 

n(7) = n(/)t (1.9) 

In particular, if / is real, then 

n(/) = n(/)t (i.io) 

in other words fi(/) is symmetric. The inverse of Weyl map, the so called 
Wigner map [2], is the map which associates to each operator a function of p 
and q, it is defined by 

Q-\F)(p,q) = — i— f e-^ a P+^Tr(FW\a,(5))dadP (1.11) 
(2nh) z J R2 

with F G B(T-L). In fact, combining (11. 7p and (11.111) we obtain: 

n-\n(f))(p, q ) = ^_^d«dM (ap+ ^ )r rr(o(/)wt (a) ^ (1 . 12 ) 

Tr(n(f)W(a,(3)) = Tr U^f^W^d^W^a, P)\ (1.13) 

The integral of the previous relation is in d£ drj, assuming that the trace and 
the W^(a, (3) can be switched with the integral and using the identity : 

TTe ± {aP+m= I e U«P+P<i) dpdq = (2nh) 2 5(a,P) (1.14) 
Jm 2 

we obtain 

Tr (Q(f)W(a, (3)) = 2nh [ f(C,vM ~ «, V ~ P)d£*n (1-15) 

Jr 2 

Therefore (11.121) becomes: 

n- 1 (n(f))(p,q) = -L [ f(a,0)e&°* + MdadP (1.16) 

The left side is exactly the inverse Fourier transform of f(a,(3). Finally we 
obtain: 

n- 1 (n(f))(p,q) = f(p,q) (1.17) 
this shows that the Wigner map is actually the inverse application Weyl map 

«nu). 

For simplicity, this construction used functions defined on IR 2 but can be 
extended to R n without difficulties. In general using an n-dimensional base is 
possible to define a one to one correspondence between a function of p and q 
with an operator using their expansion in each space 
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1.2 Moyal product 

Using Weyl map it is possible to define a new product in order to describe 
a quantum system, the Moyal product or star product. This new product is 
no longer commutative and is obtained combining Weyl map applied on two 
elements of the continuum space, with the Winger map applied on the product 
of the previous transformed elements. 
Formally the Moyal product [22] is defined by 

f* g = n- 1 (n(f)n(g)) (lis) 

This star product is associative but non- commutative and satisfies the condi- 
tion 

(af)-kg = a(f-kg) = f-k(ag) (1.19) 
with a G C. In addition, it satisfies the condition 



f*9 = 9*f, (1-20) 

the Leibniz rule 

with i G {1,2, ... ,n} and the condition 

f(x) -k g(x)dx — / f(x)g(x)dx (1.22) 
from which follows that the integral has trace property: 

f(x) * g(x)dx — / g(x) * f(x)dx (1-23) 



It is possible to show [29] that, for / and g analytic, the star product can be 
written as 

(f*9)(q,p) = /(g,p)e f V**-**J^( g , p ) (1.24) 



or equivalently 



(f*g)(q,p) = Y]f(q,p)- 



n=0 



ih I _d__d d_d_ 

2 I dq l dpi dq l dpi 



g{q,p) (1.25) 



where the arrows point the direction in which the partial derivate acts. 
The expression is well defined only on a quite small set of functions, however 
can be obtained some integral expressions with larger domains 
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(/ * g)(u) := [ [ f(v)g(w)e^> Q - 1 <^dvdw (1.26) 



n .71Bn 



Equivalent formulas are 

1 

P) 



(f*g)( u ):=—L- I I f(u + v)g{u + w)ef v - @ - x - w dvdw (1.27) 



or 

tf * 9)(v) r4s=Z f I f(u + ^e-v)g(u + w)e i ^dvdw (1.28) 

A very useful remark is that Moyal product becomes the usual one in the 
limit H — > 0, in other words Moyal product can be viewed like a deformation 
usual product in the deformation parameter fx. Therefore, we can interpret the 
quantum mechanic phase space as a deformation of classical mechanic phase 
space obtained substituting the point wise product with Moyal one. This 
product will be used in the next chapter to construct a non-commutative M@ 
vector space the so called Moyal plane. 

It is possible to give some matrix basis for the algebra of functions on a two 
dimensional phase space with the star product. Defined the functions (taking 

n = i) 

a=^ 2 {q + iq) a = j-(q - ip) (1-29) 
for a generic function results 

a*f = af+% f*a = af-% (1-30) 
and similar relations for a. The function 

2 2 

Vo = 2e^ (1.31) 
is Gaussian and has the useful property that 

^0*^0 = ^0 (1.32) 

Defining the functions 

<p m n = — / == = , , a n *<lo*a m (1.33) 
V2 n+m n!m! 

and using the relations ( 11.30]) and ( II .32]) it is easy to show that 

Therefore the yj's can be used as a base [29 j for the deformed algebra and the 
star multiplication becomes the usual multiplication of (infinite) matrices. 
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1.3 Spectral geometry 

In this section we introduce the argument of Connes' spectral geometry which 
is the non- commutative generalization of a usual geometry on a manifold. In 
particular we will introduce some basic concepts like non- commutative infinites- 
imals and the Dixmier trace as algebraic generalization of the usual infinites- 
imals and of the integral Offi l31\ \2J$ . A fundamental role will be taken by 
the generalized Dirac operator, with which we will define the metric of a non- 
commutative space and the spectral action of a Connes-Lott Model I5ffi [5lf . 

1.3.1 Non-commutative infinitesimals 

To define the Diximier trace we need some facts about compact operators [J7]. 
Let us recall that an operator A belonging C*-algebras of bounded operators 
over an Hilbert space namely A G 13(H) , is called of finite rank if the orthog- 
onal complement of its kernel is finite dimensional. Roughly speaking, even if 
the Hilbert space is infinite dimensional, such operators are finite dimensional 
matrices. Beside an operator A on TL is compact if it can be approximated in 
norm by finite rank operators. An equivalent way to characterize a compact 
operator A is that Ve > exists a subspace of finite dimension E C H which 
satisfies || A \ E ± ||< e. 

Taking in account the previous characterization of compact operators they can 
be seen in some sense small, so they are good candidates to be infinitesimals. 
Calling 1C(H) C 13(H) the subset of the compact operators, the size of the 
infinitesimal A G K(H) can be defined by the rate of decay of a sequence 
a n (A) when n — > oo. Where {a n (A)} are the non vanishing eigenvalues of the 
operator \A\ arranged with repeated multiplicity. 

The infinitesimals of order a are all the operators A G K.(H) such that it is 
possible to to construct a sequence satisfying: 

a n (A) = 0(n~ a ) as n — )• oo, 

(1.35) 

3 C < oo | a n (A) < Cn~ a ,\fn > 1. 

Remark: 

The algebra Yii^K) is the only norm closed and two-sided ideal when H is separable 
and it is essential, therefore is the largest two-sided ideal in the C*-algebra B(H). Since 
the identity operator I on an infinite dimensional Hilbert space is not compact the algebra 
is not unital. Beside the defining representation of K(H) is the only, up to equivalence, 
irreducible representation of K(H) : in fact it is Morita equivalent to the algebras of finite 
rank matrices of complex numbers. 
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1.3.2 The Dixmier Trace 

The Dixmier trace is a trace on a space of linear operators on an Hilbert 
space larger than the space of trace class operators. This trace is defined 
in order to have the infinitesimals of order 1 in its domain but the higher 
order infinitesimals have vanishing trace. With this requirements the standard 
trace, with domain in the two-sided ideal C l of trace class operators, is not 
suitable. The usual trace, defined as Tr A := J2 n (A£ n ,£ n ), for any A G C 1 , is 
independent of the orthonormal basis £ n , n G N of %. In the case of a positive 
and compact operator A, the eigenvalues becomes positive and the definition 
turn to be TrA := Yl^=o a n{A). In general, an infinitesimal of order 1 is not 
in C 1 , since we can not ensure the convergence of its characteristic values, we 
can just say that a n (A) < CrT x for some positive constant C . But we are 
sure that C 1 contains infinitesimals of order higher than 1, additionally for 
positive infinitesimals of order 1, the standard trace is at most logarithmically 
divergent since Tr A = J2n=o a n{A) < ClnN. 

The aim of the Dixmier trace [H] is to extract the coefficient of the logarithmic 
divergence, enlarging the domain of the usual trace. Let £( 1,0 °) indicate the 
ideal of infinitesimal of order 1. As a first try, if A G Z^ 1 ' 00 ) is positive, we 
could think to define a positive functional using the limit 

1 N ^ 

lim — - Va n M) (1.36) 

n=0 

But this definition is affected by two big problems: it is not linear and the 
lack of convergence. Dixmier proved [IB] that exists an infinite number of scale 
invariant linear forms lim w on the space Z°°(N) of bounded sequences. For each 
such form lim. there is a positive trace on the positive part of £( 1,0 °) defined 
by 

i tf-i 

Tr w (A) = lim.— £ a n (A), VA G A > (1.37) 

n=0 

This trace is called the Dixmier trace, it is invariant under unitary transfor- 
mations since the eigenvalues a n (A) and the sequence {7^} are invariant as 
well. It satisfy the usual trace properties and 

Tr w (A) = 0, if A is of order higher than 1 

This last statement follows from the fact that the space of all infinitesimals of 
order higher than 1 forms a two-sided ideal and its elements satisfy 



a n (A) 



n 



, or na n (T) — > as n — > 00 



1.38) 
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The higher order sequence {7^} converges to zero vanishing the Dixmier trace. 
In many practical examples in physic, like gravity and Yang-Mills theories, the 
sequence {7^} converges itself. In these cases, from the property of Dixmier 
trace, the trace is given by (11.36P and does not depends on the choice of linear 
form u, this type of operators are called measurable [28] 1^9] - 

1.3.3 Spectral triples 

We now define the spectral triple, it will be the core object of the Connes- 
models. A spectral triple (A,H,V) is defined by an involutive algebra A 
represented as an algebra of bounded operators on the Hilbert space H and 
with a self-adjoint operator T> = T>^ on "H with the properties: 

1. The resolvent {T> — A) -1 , A G f , is a compact operator on TL; 

2. [V, a]=Va-aV G B(H), for any a G A. 

The first statement tell us that the self-adjoint operator T> has a real dis- 
crete spectrum of eigenvalues A n and each eigenvalue has a finite multiplicity. 
Additionally, |A n | — > 00 as n — > 00; since (T> — A) -1 is compact, it has char- 
acteristic values a n (V — A) -1 — > 0, from which |A n | = a n (|P|) — > 00. The 
second condition can be relaxed to be satisfied only for a dense sub-algebra of 
A. A spectral triple (A,T-L,V) is said to be of dimension n > if |£>| _1 is an 
infinitesimal (in the sense of (I1.35P ) of order 1/n or, in other words, if |P| _n 
is an infinitesimal of order 1. 

The dimension of spacetime is a local property and it can be equivalently found 
from the asymptotic behavior of the spectrum of the Dirac operator for large 
eigenvalues. Ordering the eigenvalues,- • • A„i < A n _i < A„, + i ■ • • , the Weyls 
spectral theorem states that the eigenvalues grow asymptotically as ri 1//dimM . 
The local property of spacetime are encoded in the high energy part of the 
spectrum. This agreement with our intuition from quantum mechanics and 
motivates the name spectral triple. In following pages many regularity condi- 
tions on elements of A will be defined [31 J using only the operator V and its 
modulus I'D I and T> will be a generalization of the usual Dirac operator on an 
ordinary spin manifold, for simplicity we will call it just Dirac operator. 
In order to find the analogue of the measure integral, the Dirac operator will 
play a role for definition of the volume. With such notion we can define the 
integral of any b G A, for a n-dimensional spectral triple, by the following 
formula: 



Where the operator \T>\ n is just used to bring the bounded operator b into 
£(i,oo) ) in this 

way the Dixmier trace is well defined. In the previous definition 




(1.39) 
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operator the |"D|~ n is the analogue of the volume form of the space and can be 
proved that the integral (11.391) is a non-negative normalized trace on A given 
a general spectral triple (AyHyV) |3T] . 

Now let us introduce a more structured spectral triple: the even spectral 
triple. The object composed by five items (A, H, T>, J, x) form what Connes 
calls |35j an even real spectral triple. Again A is a real associative involution 
algebra with unit, represented faithfully by bounded operators on the Hilbert 
space "H, V is an unbounded self adjoint operator on "H. In addition, J is an 
anti-unitary operator and x a unitary one. They fulfill the following properties: 

1. J 2 = — 1 in four dimensions (J 2 = 1 in zero dimensions) 

2. [p^^Vp^V- 1 ] = for all a, a G A 

3. DJ = JD 

4. [V,p(a))] is bounded for all a G A and [[£>, p(a)} , Vp{a)V- 1 } = for all 
a, a G A 

5. x 2 — 1 an d [Xj p( a )] = for all a G A 
6- Jx = XJ 

7. V X = -XV 

Furthermore, can be required some other properties: Poincare duality, regular- 
ity, which grands that our functions a G A are different iable and orientability, 
which connects the volume form to the chirality. The fourth property is called 
first order condition because it grants that the Dirac operator is a first order 
differential operator and property 5 allows the decomposition H = T-Ll® Hr. 
These properties were promoted to axioms by Connes defining an even real 
spectral triple justified by his Reconstruction Theorem [36]. Considering an 
even spectral triple (AyT-LyV, J,x), with commutative algebra A, then exists 
a compact Riemannian spin manifold M of even dimensions whose spectral 
triple (C 00 (M),L 2 (S),^,C,7 5 ) coincides exactly with (A,n,V, J,x). Where 
C°°(M) is the algebra of the infinity derivable functions on M, L 2 (S) is the 
space of square integrable functions on the spinor space S, $ = ydi is the 
standard Dirac operator, C is the charge conjugation operator and 75 is the 
usual chirality operator 75 = 71727374 • • -Jn- 

In this theorem are contained a lot of informations about the role of the 
Dirac operator. Beside, describing the dynamics of the spinor field, the di- 
mension of spacetime and its integration, the Dirac operator $ encodes its 
Riemannian metric and its differential forms. The metric can be reconstructed 
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from the spectral triple by Connes distance formula; a point x G M is re- 
constructed as the pure state and the general definition of a pure state of 
course does not use the commutativity. A state 6 of the algebra A is a lin- 
ear form on A, that is normalized 5(1) = 1 and positive 5(a*a) > for all 
a G A. A state is called pure if it cannot be written as a linear combination 
of two states. In the commutative case, there is a bijective correspondence 
between points x G M and pure states 5x defined by the Dirac distribution, 
S x (a) := a(x) = J M S x (y)a(y)d 4 y. The geodesic distance between two points x 
and y is reconstructed by the Connes distance formula: 

sup{|k(a) - 5 y (a)\; a G C°°(M) | || [V,p(a)] ||< 1} (1.40) 

For a general even spectral triple the operator norm || [D, p(a)\ ||, in the dis- 
tance formula, is bounded by axiom. 

For example consider the circle M = S 1 of circumference 2n with Dirac op- 
erator d = id/dx. A function a G C°°(S 1 ) is represented faithfully on the 
space of square integrable functions ip G L 2 (S' 1 ) by pointwise multiplication, 
{p(a)t/))(x) = a(x)ip(x). The commutator [d, p(a)] = ip(a!) is bounded and we 
have already seen it in quantum mechanic. The operator norm in this case is 

II [d,p(a)} ||:= 2tt sup^d, p(a)]ip\/\ip\ = supja'(a;)| (1.41) 

Where \ijj\ 2 = ip(x)ip(x)dx. Using (11.401) we recover the standard distance 
between two points x and y on the circle: 

sup{|a(:r) — a(y)\ \ sup x \a'(x)\ < 1} = \x — y\. (1-42) 

It is important to note that Connes distance formula works for non-connected 
manifolds too, even for discrete spaces of dimension zero or collection of points. 
The last ingredient that we need are the differential forms and again they can 
be recovered using the Dirac operator by an analogy with quantum mechanic. 
The differential form of degree one, like da, for a function a G A are recon- 
structed as —i[D,p(a)]. For example a 1+1 dimensional spacetime, da is just 
the time derivative of the observable a and is associated with the commutator 
of the Hamilton operator with a. Higher degree differential forms are obtained 
by multiple application of the commutator with T>. 

We define a non-commutative geometry by a real spectral triple, which con- 
tains all the geometric informations, with non-commutative algebra A . 

1.3.4 Spectral Action 

[4*0] The axioms of the spectral triple allow us a change of point of view. A 
quite suited analogy is the Fourier transform, in which the points x G M of 
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a geometric space are replaced by the spectrum S C 1 of the operator T>. 
In fact, we can forget about the algebra A in the spectral triple (A : H : V) 
and focus ourself only on the operators V, 7 and J acting in H and we are 
able to characterize this data by the spectrum E of T> which, for the even case 
E = — E, is a discrete subset with multiplicity of R. So we can argue that all the 
physical informations, therefore the physical action, only depends the spectrum 
of V. The next step is to look for the existence of Riemannian manifolds 
which have the same E, namely isospectral and in general not isometric. The 
previous hypothesis is much stronger than the invariance of the action under 
diffeomorphisms of general relativity. This approach has the virtue of not 
require the commutativity of the algebra A in order to apply this principle 
to a physical action. Indeed, the spectral triple axioms require just the much 
weaker condition between the algebra and the opposite algebra: 

[a, b°] = Va, b G A with b° = Jb* T 1 (1.43) 

Analyzing, in the Riemannian case A = C°°(M), it is possible to construct 
an isomorphism between the group of diffeomorphism of the manifold Diff(M) 
and the group automorphisms of the algebra Aut(^4). To each diffeomorphism 
G Diff(M) it is associated the algebra preserving map : A — > A given by: 

Mf) = f°<t>~ 1 VfeA = C°°(M) (1.44) 

This association is in general true, the group Aut(*4) of automorphisms of the 
algebra A is the generalization of the diffeomorphisms to the non- commutative 
geometry (A, %,"£>). It is important to notice that in the general case there is 
a non trivial normal subgroup of the group Aut(^4) 

Int(^) C Aut(^) (1.45) 

where Int(^4) is the group of inner automorphisms; n is called inner if exists a 
unitary operator u G A such that 7r(a) = uau* Va G A. 

The subgroup Int(>4) of inner automorphisms is a normal subgroup and the 
group Aut(^4) of diffeomorphisms falls in equivalence classes under Int(^4). 
This induces a natural foliation into equivalence classes in the space of metrics. 
The internal fluctuations of a given metric are given by the formula, 

V = V + A + J A J' 1 , Ae&iA) (1.46) 

where 

tt\A) := {A = E^Do, k], I a h k G A and A = A*} (1.47) 

Applying the previous formula to (A,'H,'D ), where V is the un-fluctuated 
V, the fluctuations does not affect the representation of A in V. but perturbs 
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the operator T> by ( I1.46P where A is an arbitrary self-adjoint operator in % of 
the form A = £aj[X>o,&i] aj,6j G A. The fluctuated Dirac operator continues 
to satisfy the axioms and in the commutative case (i.e. Riemann) the group of 
inner automorphisms Int(^4) = {1} is trivial, as a consequence the fluctuations 
are trivial too T> = T>q 

The action of hit (A) (where the tilde stands for taking into account the action 
of automorphisms on the Hilbert space ~H) on the space of metrics is restricted 
on the above equivalence classes and is given by: 

u£u* = uu*°£, A u[V, u*) + uAu* (1.48) 



From the previous described properties of a general real spectral triple follows 
that it can be used to define a gauge theory. The gauge fields are recovered 
from the inner fluctuations of the Dirac operator and the gauge group is given 
by the unitary elements in the algebra, we still need to define an physical 
action. Let (^4, H, V, J) be a real spectral triple, given the fluctuated operator 
V, a positive even function \ an d a cut-off scale A, it is possible to define a 
gauge invariant spectral action for the bosonic sector: 

S b [A\:=Tr X (j^ (1-49) 

The cut-off parameter A is used to obtain an asymptotic series for the spectral 
action, in this way the physically relevant terms will appear as the coefficients 
of an expansion in positive power of A. For the fermionic sector, we can define 
a fermionic action in terms of ip £ H and A G fi 1 (^4): 

S f [A,1>]:=(ii>,1Xi>) (1-50) 



In the last part of the chapter we sketch a physical application of the 
spectral action, treating the case of standard model plus gravity described by 
the action functional 

S = S E + S SM (1.51) 

where Se = j^g J R^/g d 4 x is the Einstein action and Ssm is the standard 
model action. It involves, beside the metric, other fields: bosons of spin such 
as the Higgs, bosons of spin 1 like 7, W ± and Z, the eight gluons, fermions, 
quarks and leptons. The two parts of the action have a priori a very different 
origin; the interaction of Ssm, which is governed by a gauge invariance group, 
is a priori very different from the interaction of the Einstein action which is 
governed by invariance under the diffeomorphism group. Formally the natural 
group of invariance of the functional fll.5ip is the semidirect product, 

G = U x Diff(M) (1.52) 
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Where U = C°°(M, U(l) x SU{2) x SU{3)) is the group of local gauge trans- 
formations. 

It is very useful to compare the exact sequence of endomorphisms groups of 
A, 

1 ->■ Int(.A) ->■ Aut(A) -> Out (.A) ->■ 1 (1.53) 

with the exact sequence which describes the structure of the symmetry group 
G of the action functional (11.511) . 

1 -> U -> G ->■ Diff(M) -> 1 (1.54) 

The (11.531) and (11.541) look very similar and a natural question arises: to find 
an algebra A which satisfy Aut(*4) = G condition. Taking into account the 
action of inner automorphisms of A in % given by 

£ u{u*f£ = u£u* (1.55) 

this algebra turns out to be: 

A = C°°(M) £g> Af (1-56) 

With 

A F = C ®H® M 3 (C) (1.57) 

The algebra Af is finite dimensional and H C M2(C) is the algebra of quater- 
nions, 

/ a \ 



a (3 
\-0 a j 



a,(3eC (1.58) 



Next to the algebra A, we need the action of A in H and a suitable Dirac 
operator operator T>. The algebra A (I1.57P is a tensor product of two algebras 
which corresponds to a product of spectral triples given by: 

(AH,P) = (C°°(M),L 2 (M,5), $ M )®(A F ,U F ,V F ) (1.59) 

with 

H = L 2 (M, S) ®U F i V = $ m ®1 + ^®V F (1.60) 

The algebra Af is finite dimensional so the dimension of the corresponding 
space is and T-L F must be finite dimensional [51]. The elementary fermions 
provide a natural candidate for % F and a finite Hilbert base can be labeled by 
elementary leptons and quarks. For instance, for the first generation of leptons 
we get 6l, e#, u L , e L , e#, v L as the corresponding basis. The helicity operator 
■jf is given by the usual 75 and distinguishes the left handed particles and 
right handed ones. For quarks one has an additional color index, g,r,b. The 
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real structure J is just such that J f = f for any / in the basis. Additionally, 
the algebra Af has a natural representation in Tip and: 

ab° = b°a Va, b G A F , b° = JUT 1 . (1.61) 

The finite Dirac operator acting in the finite dimensional Hilbert space Af 
which fulfill the spectral triple axioms fl 1 . 43 1) is: 



\° 9 1 



(1.62) 



where Y is the Yukawa coupling matrix. 

Now we are able to determine the internal fluctuations using (ll.46p . computing 
the internal fluctuations of the above product geometry M x F, skipping the 
calculus for brevity, we recover the bosons [58, 50]. In fact, the internal fluctu- 
ations are parametrized exactly by the bosons 7, W and Z the eight gluons 
and the Higgs field H of the standard model. The last step is to compute 
the spectral action for the fluctuated V, the calculus are quite cumbersome 
but it is possible to prove [5H1 EH ED] that for any smooth cut-off function x, 
x(A) = 1 for |A| < 1, we have: 

Tr X (t) =S e + S g + S gh + S h + S c + 0(A-°°) 



The last missing contributions, the femionic sector, in terms of the operator 
T> alone are given by the equality: 



J M 



In [H] this construction was improved to include massive neutrinos and to 
solve some technical issues jl2] at the same time. 

For completeness, we end this chapter with some mathematical tools useful to 
compute (in particular conditions) the spectral action using the heat kernel 
expansions and Seeley-De Witt coefficients [HO]. For a vector bundle V on a 
compact Riemannian manifold (M, g) and a second-order elliptic differential 
operator P : C°°(V) C°°(V) of the form 

P = -{sTdpdv + + L) (1.63) 

with K^,L G r(End(V)), it is possible to find a unique connection V and an 
endomorphism E on V satisfying: 

p = VV*-£. (1.64) 
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Or locally V M = <9 M + u', where 



^ = ^K» + g^T; a ) (1.65) 

Using this u'^ and L we find E G r(End(V^)) and we can compute the curvature 
Sl^v of V: 

E := L - sTQuK) ~ iTWu + (1-66) 
:= ^(^) - ^H) - K,^] (1-67) 

Now it is convenient to make an asymptotic expansion (as t — Y 0) of the trace 
of the operator e~* p in powers of t: 

Tre~ tp ~J2t {n - m)/2 a n (P), a n {P) := / a„(x, P)^d m x (1.68) 

n>0 7M 

The coefficients a n (x, P) are called the Seeley-DeWitt coefficients and m is the 
dimension of M. Can be proved \§U\ that a n (x,P) =0 for n odd and the first 
three even coefficients are given by 



a (x,P) = (4vr)- m/2 Tr(id) (1.69) 
a 2 (x,P) = (47r)- m/2 Tr(-i?/6 id + E) (1.70) 

a 4 {x,P) = (4n)- m / 2 ^Tr(-l2R^ + 5R 2 -2R^ + 

2R fXVfXT W lVf " T - QORE + 180E 2 + 60Ef + 30^^) (1.71) 



where Rf 1 := V M V M -R and E? := V^V^-E. Considering manifolds without 
boundary, the terms E? , vanish due to the Stokes Theorem. This ex- 
pansion is very useful in some computations of the spectral action. Taken a 
fluctuated Dirac operator T> for which T> 2 can be written as (11.631) on some 
vector bundle Vona compact Riemannian manifold M and writing \ as a 
Laplace transform, we obtain 

X (V A /A) = [ ~g(t)e- W - /A2 dt. (1.72) 
Jt>o 

Using (11.681) we find that for a four- dimensional manifold the relevant terms 
of the expansion are 

Tr x(2VA) = 2k\ i a Q {V 2 A ) + 2A 2 X2 a 2 (P 2 ) + a 4 (X? 2 ) X (0) + O^ 2 ), (1.73) 

where Xk are the moments of the function x '■ 

POO 

Xk ■= / x(w)w k ~ l dw; (k > 0). (1.74) 
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8-dim spectral action 



In this chapter will be computed a spectral action starting from a non-com- 
mutative spectral triple. The feature of this particular triple is the choice 
of a 4-dimension Harmonic Dirac operator. The idea behind this construc- 
tion [45J is to relate the Dirac operator, which plays a fundamental role in 
a spectral triple, with the oscillator Hamiltonian operator. Roughly speaking, 
we look at the Dirac operator as a generalization of the Laplace operator so 
we have T> 2 rj H. Considering the spectrum u(2n + 1) n G iV of the one- 
dimensional harmonic oscillator Hamiltonian H, can be deduced that H" 1 is a 
non-commutative infinitesimal of order one. According to the previous discus- 
sion the non-commutative dimension of a spectral triple, equipped with the 4D 
harmonic Dirac operator T>\ w H4, is fixed by the order of the inverse operator 
T? -1 which is eight not four. Due to the fact that the spectral dimension is 
defined by the Dirac operator, it is connected to the phase space dimension 
and not on the one of the configuration space [54J . In order to construct such 
harmonic Dirac operator and the spectral triple we will work in the frame- 
work of the generalized n-dimensional harmonic operator. Therefore, will be 
studied the four- dimensional case in order to construct the non- commutative 
spectral triple which is starting point for the field theory we are interested 
in. The first part of the chapter will be devoted to introduce the harmonic 
Dirac operator starting from the one-dimensional case and then generalizing 
it using the annihilations and creations operator for the bosonic and fermionic 
sectors. Having the 4-dimensional harmonic Dirac operator with harmonic 
oscillator spectrum, to implement the Higgs mechanism we will consider the 
tensor product of the non-commutative triple with a finite Connes-Lott type 
spectral triple [37] • We will fluctuate the total Dirac operator following the 
previous described standard machinery [4"0~1 of non- commutative geometry 
to get " Gauged" Dirac operator. Thus we will proceed to compute the spectral 
action in which are present two U(l)-Moyal Yang-Mills fields unified with a 
complex Higgs field. 
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2.1 Harmonic Dirac operators 

To introduce this subject we start to show the simple case of one dimensional 
harmonic oscillator, then will be introduced the general n-dimensional case. 
In one dimension the Hamiltonian of an harmonic oscillator is well known and 
is defined as 

Hi = -^- + u 2 x 2 (2.1) 
ax z 

Using the Hermite polynomials is possible to construct an orthonormal base of 
eigenf unctions of Hi in the Hilbert space L 2 (R) of square integrable functions 
on R 

/9" 



H l( p n = cu(2n + l)i; n , ^ n {x) = (f) * *^gEl e ~¥, n EN (2.2) 



where H n (z) are the Hermite polynomials. From the behavior for n — > oo 
of the eigenvalues of Hi we can infer that the inverse operator 1 is a first 
order non commutative infinitesimal. Reminding previous chapter, it means 
that after arranging the eigenvalues in deceasing order taking in account the 
multiplicities, the order of the generic eigenvalues a n is 0{n^ 1 ). From this 
evidence we can deduce that taking a Dirac operator to be V 2 m Hi, the 
non-commutative order of this Dirac operator will be 2. In order to define 
such Dirac operator we need a Clifford algebra of dimension two which can be 
represented by the Pauli matrices. In this setting a choice of a two dimensions 
harmonic Dirac operator can be: 



• d 
T>i = i<j\- — Yviwx 

ax 



\ * (i - H / 



\ 



V2ua 
2uja) 

(2.3) 

Where a = -k=(4- + ux) and = -k={4 — ux) are the usual annihilation 



ad creation operators which satisfy the commutation relation [a,^] = 1. To 
complete a spectral triple data we need to define the Hilbert space on which 
this operator acts, an algebra and its representation on H. It easy to see that 
the Hilbert space is just Hi = L 2 (R) cg> C 2 , about the algebra the simplest 
choice is the commutative algebra of the Schwartz class functions A = «S(R). 
The representation can be defined as the pointwise diagonal multiplication 
7r (/)(V'i; ^2) = (/^i,/^) with ipi,ip2 ^ L 2 (IR). This spectral triple can be 
turned in a even spectral triple finding a suitable real structure Ji, using the 
axioms of even spectral triple J 2 = e, J{D = e'VJi, JiX = e"xJ\ and the anti 
linearity, Ji turns out to be (Jii/;)(x) = a^ij){x) with the choice e = 1, e' = 1, 
e" = 1. Can be proved [57] that all the axioms of a even spectral triple are 
fulfilled, in particular the constraints about the opposite algebra are satisfied 
by this choice of J x and the operator [D, 7r(/)] = ioi7r(/') is bounded. 
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The previous construction can be generalized in n-dimensions using the 
Clifford algebra of IR 2 ™ represented on the Hilbert space C 2 ™ , in order to do it 
is very useful to consider b, the fermionic annihilation and creation operators 
with the usual anti-commutation rules: 

{6,6t} = i 

{b,b} = (2.4) 
{&t,&t} = 

With these operators we can redefine the previous one dimensional spectral 
triple as: 

£>! = -i\f2uoa ] ®b + iV2ua <g> 6 t = i^- ® (b + b ] ) + iux ® (b — b ] ) 

ax 

Hi = L 2 (M)®(C|0)©C6 t |0)) xi = 1 ® (& f 6 - 6& f ) ( 2 -5) 

Where C|0) is the space generated by the vacuum state defined as b|0) = 0. 
Now the generalization to n-dimensional harmonic oscillator becomes straight- 
forward if we considering n-dimensional fermionic annihilation and creation 
operators b^, b\, and n-dimensional bosonic annihilation and creation opera- 
tors £fy, a J, satisfying for /jl, v — 1, • • • , n: 

[a-n, a v ] = [a*, aj,] = 0, [a M , a£] = 5^ (2.6) 
{&„, 6,} = {fet , fet } = o, &t } = ^ (2.7) 

Where a M = t/^O-^ju + 9^), = -75^(^2^ — d^). The generalization of the 
Dirac operators (I2.3P using this operator is: 

V n = -iV^jS^al ®K + iV^uS^a^ ® b\ = ® (&„ + &t ) + zwa^ ® (&„ - b\) 

(2.8) 

summed over repeated index. In analogy with the one dimensional case we 
can define the Hilbert space on which the Dirac operator (12. 8p acts starting 
from the vacuum state by subsequent applications of the fermionic creation 
operators bj, on the vacuum, using the anti-commutation relations (12.71) and 
the definition of vacuum state &|0) = 0. We call this space A(C n ) and therefore 
the Hilbert space is l-i n = L 2 (R") <g> A(C n ). Beside, we can define a grading 
operator Xn as: 

n 

Xn = I®IlM- 6 4) ( 2 - 9 ) 

n=i 

An interesting feature of the Dirac operator (12.81) is the preserving of the sums 
of excitations of bosons and fermions in the Fock space, in other words is 
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super-symmetric. However, the super-symmetry will no longer holds for the 
total spectral triple due to the non-preserving behavior of algebra that will be 
chosen. Using the relations ( 12. 6ft - ( 1277)) we can compute the square the Dirac 
operator ( 12. 8p as: 

V 2 n = 2ua\a? ® I - 2wl ® b\W = 2coN B <8> I - 2wl <g> N F (2.10) 

Where Np and Nb are the number operators. In this form it easy to see that 
P 2 , being a "difference" between fermionic and bosonic number operator, has 
only one zero mode corresponding to the vacuum state. For practical reasons 
it is convenient write D\ as: 

Vl = {a\ a v + a u aD® I- 2wl ^ (&*&„ - b v b\ ) = H n ® I +u;®£ n (2.11) 

where in H n we can recognize the harmonic oscillator Hamiltonian and the 
spin operator £„. The universality property of the Clifford algebra grants 
the existence of an isomorphism between the 2n- dimensional Clifford algebra 
and the Hilbert space H n = L 2 (R") <g> C 2 ™. In this representation the Dirac 
operator is: 

V n = zTO M + wI^ + % (2. 12) 

Where P« turns to be 1^ = (b^ + &t), I^ +n = i(&t - 6 M ) which satisfy the 
relations: 

T a r b + T 6 r a = 25 ab with a, & = 1, • • • , 2n (2.13) 
Beside, the grading operator is represented as: 

xn = H) n (-i)^ g» r\ • • • r 2n (2.14) 



2.2 An harmonic spectral triple for the Moyal 
plane 

In the framework of non-commutative field theories on 4-dimensional Moyal 
plane has been proved [HI H3j that the introduction of an harmonic oscillator 
term makes a 4 -model on 4-dimensional Moyal plane renormalizable. Such 
oscillator term can be written as: 

H m ^--^- + n%^ + m 2 (2.15) 
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where x 
9 = i6a with a 



/. • 2(0 ) [ivX 



can be chosen as two copies of the Pauli matrix 
<r 2 ® I2 or explicitly: 



/ 



e 



0000 
-#000 
ooo# 



\ 



V 







6 e 



(2.16) 



-0 0/ 



With this choice we have _1 = — Quantum mechanics tell us that in the 
Hilbert space L 2 (IR 4 ) exists an ortho normal basis if) 8 , s G M 4 of eigenfunctions 
of H m with eigenvalues 



A s (m) 



4Q 
T 



6m 2 



s G N 



(2.17) 



The inverse extends to a selfadjoint compact operator on L 2 (R 4 ) with 
eigenvalues \~ 1 {m). If we look at the trace the operator H~ 4 we find: 



Tr(F-) = X> + 3)(n + 2)(n + l)(A„(m))< 



(2.18) 



n=0 



Which is derived simply from the number of possibilities to express s as a 
sum of four ordered natural numbers. This means that the trace exists or in 
other words H~ 4 belongs to the Dixmier trace ideal L( 1,00 )(L 2 (IR 4 )) of compact 
operators. Nevertheless, the main object of non-commutative geometry, which 
determine the spectral dimension, is the Dirac operator. The relation T> 2 = H, 
which essentially states that each degree of freedom in the configuration space 
contributes to the spectral dimension twice, implies that the 4-dimensional 
Moyal space has spectral dimension 8. 

From the previous section, we can define a proper Dirac operator just 
considering the 4-dimensional case obtaining a Dirac operator built from a 
8-dimensional Cifford algebra: 



(2.19) 



Here, the Tk G Mi 6 (C), k = 1, 8 are the generators of the 8-dimensional real 
Clifford algebra, satisfying 



(2.20) 

1 ® C 16 of square integrable 
spinors over 4-dimensional euclidean space. Accordingly with (12.131) for ip G 
%4 we obtain: 



TaT^ + TiT k = 25 k il 
We take the Hilbert space H A = L 2 (M 4 , S) = L 2 (M 4 ) 



V\i) = ((-A + ft%^)I + E) ip , E := -itt(Q- l y v [T ^ V 



(2.21) 
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with A = <9 M <9 M . As algebra we chose the Moyal algebra M@: 

A 4 = R% = (S(K 4 ),*) (2.22) 

where (iS(R 4 ), *) is the algebra of the Schwartz functions on IR 4 , with the Moyal 
product 

(/*»)(*) = / / dS-0yJ{x+ l -Q-k)g{x + y)e^ , f,geA 4 (2.23) 

The representation of the algebra A± on % 4 is by component-wise diagonal 
Moyal product [55] * : A± x — >■ H4. The Moyal product can be extended 
to constant functions using another representation of the product with the 
integral representation of the Dirac distribution. Taking in account, for smooth 
spinors, the identity 2x^"ip = x-kip + ip-kx and the relation 

[x v ,f] it = ie v % (2.24) 

we compute the commutator of that action with the Dirac operator 

z>4(/*V0-/*(z>4V0 



= i(r^ + fir^ +4 )(a At /)^^- 



(2.25) 



Just the four-dimensional differential of / appears, no x-multiplication, this 
differential is represented on by 7r(<ia; M ) = + fir^" 1-4 . Furthermore the 
previous commutator confirms that (A^l-L^V^) satisfy the mak@ axioms of 
spectral triple, in fact the commutator is bounded and due to its commutation 
with Moyal multiplication, order-one condition is fulfilled. 

Now we introduce a very useful relation connected to the heat kernel type 
expansion associated to a regular spectral triple [57J. This relation will be 
used later in order to compute the spectral action. Considering a regular 
non-unital spectral triple (A,T-L,T>, J) and two pseudo-differential operator 
A G ^o(A) A\ G ^i(A) of order respectively and 1. we consider the 
following decomposition: 

4 

e -t{vi + A 0+Al ) = J^i-iyEjit) - t 5 R (2.26) 

j=0 



^dn the appendix A a unitalised Moyal algebra is introduced in the frame of non compact 
spectral triples 

2 Orientability axiom and Poincare duality will be not considered 
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Using Duhamel principle [56] 

e -t(A+B) = e -tA _ t f 1 dg e -st(A + B) Be -(l-s)tA 
JO 



we can identify: 

E (t) = e- w2 
EM : 



J a, 

(2.28) 

and 

R= f d Sl ds 2 ds 3 ds 4 ds 5 e- Slt{v2+Ao+Al) (A + Ai)e^ 82 ~ 8l)wi ''{A + A,) 
Ja s 

xe-^- S2)tv2 (A + AJe'^-^iAo + A^e^-^ (A + A^ 1 ' 8 ^ 

(2.29) 

The domains of the integrals Aj are the j-simplex: 

j 

Aj = {s e R j ; < si < s 2 < si < ■ ■ ■ < Sj < 1} ~ {s e R j+1 ; s< > 0, ^ s { = 1} 

i=0 

(2.30) 

Taking in account the relation: 

e~ w \A] = j' ds^ {e~ tsV2 Ae- l ^ v2 ) = -t j" dse~ tsV2 [V 2 ,A] e -*Ci-P a 

(2.31) 
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we can rearrange the Ej collecting the heat operators as follow: 
E 2 (t) 

= [ e-^ tv \A + A 1 )e-^-^ tv2 (A + A 1 )e-^ tv2 ds 1 ds 2 
Ja 2 

= [ e~^ w \A + A 1 fe^ 1 -^ tv2 ds l ds 2 -t [ f dr(s 2 - s^e^ 2 

Ja 2 Ja 2 Jo 

x(A + A 1 )e-^- s ^ rtv2 [V 2 ,A + A t ] e -(s 2 - si )(i-r)t^ ' e -(i-«)tD= 

= / e-^ w2 {(A + A 1 ) 2 -t( S2 - Sl )(A + A 1 )[V 2 ,A + A 1 ]} 
Ja 2 

xe- {1 - Sl)w2 d Sl ds 2 + t 2 [ f dndr 2 ri(s 2 - stfe' 81 ™ 2 (A + A x ) 

J A 2 JO 

= [ e -i*» a {(A + A 1 ) 2 -t(s 2 -s 1 )(A + A 1 )[V?,A + A 1 \ 

J A 2 

+ t ^{s 2 -s l ) 2 {A, + A l ) [V 2 , [V 2 ,A + A 1 ]] }e-^ tv2 d Sl ds 2 
-t 3 [ f d Sl ds 2 r 2 ir2 (s2 - si) V^Vo + AJe-^-^^™ 2 

J A 2 Jo 

x [V 2 , [V 2 , [V 2 ,A + A,]]] e-^-^-^^e-^-^ds^ (2.32) 

The last integral is a bounded operator which tents to zero for t — > 0. 
Applying the same procedure for E 3 we get: 

E 3 (t) 

= [ e^ w2 {(A + A 1 ) 3 -t(s 3 -s 2 )(A + A 1 ) 2 [V 2 ,A + A 1 ] 
Ja 3 

-t(s 2 -s 1 )(A + A 1 ) [V 2 ,A + A 1 ] (A + A 1 )}e-( 1 - S ^ tv2 ds 1 ds 2 ds 3 

+t 2 [ [\s 3 - Sl )V^Vo + Al ) 2 e-^~ s ^ w2 
Ja 3 Jo 

x [V 2 , [V 2 ,A + A 1 \] e-^-^-^^e-^-'^dsids^Ssdnd^drs 
+t 2 [ [\s 2 -s 1 )(s 3 -s 1 )e- s ^ 2 (A + A 1 )[V 2 ,A + A 1 ] 

J A$ Jo 

Xe -( S3 - Sl )^ 2 [p2 > Aq + A J e -(^- Sl )(l-r)^ e -(1^3)^ 2 rfsirfs2rfs3rfr 

xe -( S2 - Sl )(i-r 1 r 2 )^ e -( S3 - S2 )^ 2 (Ao + AOe-^ 1 - 83 ^ 2 ^^^^ (2.33) 

Even in this case, the integrals multiplied by t 2 give rise to a vanishing bounded 
operator for t — > 0. Now using the cyclic property of the trace we can compute 
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the traces of E 2 over A 2 and of E 3 over A 3 , which considering only the leading 
terms are: 



Tr(t 2 £ 2 (t)) 

= tr({^(A + A) 2 - ^(A + A) [V 2 , A + Aj 

+ ^(A) + AO [P 2 [V 2 , A + AJ] )e-^ 2 } + O(Vt) 

= tr({^(A + Ax) 2 - ^(A [V 2 ,A 1 ]+A 1 [V 2 ,A ] + A, [V\A,]) 

+ t ^A 1 [V 2 [V 2 ,A ]} )e^ 2 }+0(V~t) (2.34) 



and 



Tr(t 3 E 3 (t)) 
t 3 



= tr({ - -(A + A) 3 + -(A + A) 2 [V 2 , A + Aj 
+^(A + A) [V 2 , A + AJ (A + A 1 ))e-*° 2 } + 0(Vf) 
= tr({ - |(A A 2 + AxAqAi + A 2 A + A?) 
+^A 2 [V 2 , Ai] + ^ [X> 2 , A x ] A^e-^ 2 } + (2.35) 
In the exactly same way we obtain for £4: 

Tr(t 4 £ 4 (t)) = t T (^(A + A) 4 e- w2 ^+O(Vt) 

= tr (^e-^) + O(Vt) (2.36) 

Summing all contributions together we finally arrive at: 

Tr(e~ <(I)2+A(,+Al) ) 

= Tr({l - t(A + A) + ^(A 2 + AxA + A A, + Aj) 

-|(A [I? 2 , Ax] + Ax [D 2 , A ] + Ax [V 2 , A x ] + A A 2 + AxA Ax + AjA 

+A 3 ) + ^(Ax [D 2 [P 2 , Ao]] + 2A\ [V 2 , A,] + Ax [V 2 , A,] A, + A 4 )e^ 2 } 
+0(Vt) (2.37) 
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2.3 8-dimensional Higgs model 



Following the Connes-Lott models, in order to implement the Higgs mech- 
anism, we consider the total spectral triple as the tensor product of the 8- 
dimensional spectral triple (A4, H4, P4, Tg) with the two point Connes-Lott 
like spectral triple (C <g> C, C 2 , Mai). The total Dirac operator of the product 
triple is: 

V T = £> 4 ® I + T 9 ® Ma x (2.38) 

Or explicitly: 

( „ ^\ 



Vn 



V A MT g 

MTg V 4 



(2.39) 



The algebra becomes At = A$ © A± and acts by diagonal star multiplication 
(I2.23P on T-Lt = Ha ® W&- The fluctuated Dirac operator is found using T>a = 
T>t + Yjidi\DT,bi\ with a^bi G At of the form (f,g), the computation of the 
commutator T>t with (/, g) gives: 



p T ,(f,g)} 



z(r^ + nr^)L*(dj) Mr 9 L*(/ - g ) 
y MT 9 L*{g-f) i(r» + nr^)L*(d,g) J 



(2.40) 



L*(f)ip = f * ip is the left Moyal multiplication. From the commutator we 
deduce that the form of selfadjoint fluctuated Dirac has to be: 



d 4 + (r„ + qt^l^) r 9 L*{ip) 

t 9 l^) v 4 + (r M + nr^L^B") 



\ 



J 



(2.41) 



Where G A± is the Higgs complex field and A^, G A4 are real fields. The 
spectral action computation needs the square of T>a- 



T> 2 



f {Hi + l,(0 * <p))i + s + f a z(r M + nr^TgL^D^) ^ 



with 
F A 



y i(r ^nr^r 9 L*(D»<f>) (h 2 + u(<p*<p))\ + ^ + f b j 

(2.42) 

d^-iAp + t + if + Bp (2.43) 

{v A , (r„ + vr^L^A")} + (r„ + fir M+4 )(r, + ^r„ +4 )L*(^ * a v ) 

{L^A"), id. + n 2 M.{x^)} + (1 + n 2 )L,(Af, * ^) 



+ % 



Q^, r„] + ^ 2 [r M+4 , r„ +4 ] + nv,v v+ A l*(*T), (2.44) 
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(M.(x M )^)(a;) = x^ip{x) is ordinary pointwise multiplication and F B is ob- 
tained just replacing A with B. We can recognize in previous expression the 
field strength = d il A v - d v A^ - i(A fM *A V -A V * 



2.3.1 Spectral action 

Recalling the spectral action principle, the bosonic action can be defined ex- 
clusively by the spectrum of the Dirac operator. The general form for such 
bosonic action is: 

S{V A ) = Tr X (^) (2.45) 

Where \ is a regularization function \ : R + — > R + for which trace exists. 
The trace in (g35D is defined on B(L 2 (R 4 )) by 

Tr(A) = dx A(x,x) (2.46) 
Jr* 

together with the matrix trace including the Clifford algebra. By Laplace 
transformation one has 

/•oo 

S{V A )= / dtTr(e- w *)m (2.47) 
Jo 

where x is the inverse Laplace transform of xi 3 )-, 

POO 

X {s) = / dte~ st x(t)- (2-48) 
Jo 

The trace in (I2.47P is given by: 

Tr(e- W *)= f dxtv(e- w ^)(x,x) (2.49) 
Jr* 

Assuming the trace of the heat kernel e~ w ^ has an asymptotic expansion 

oo 

Tr{e- W *) = a n(V 2 A )t n , 5 G N (2.50) 

n=— 5 

we obtain replacing the previous expansion into (12.471) 

OO /»OG 

S{V A ) = V a n {V\) / dt t n x(t) (2.51) 
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To compute the integrals we have to consider separately the cases n N and 
n£N. For n e N we have 



POO 

lim / dte- st t n x(t) = lim(-l) 



<9s n 



dte- M x(t) 



lim(-l)"^( S ) 



;-i) n x (n) (o) 



(2.52) 



For n ^ N we have 

POO 

/ ds s- n ~ l X (s) 
Jo 

In summary, 

Xn ~ 



OO />OC /*oo 

ds / dt e-^a-"- 1 ^^) = T(-n) / dt t"x(t) 
ii Jo Jo 

(2.53) 



r °°d ss -"- \(s) forn^N 
(-l)"-* x (™)(0) for n G N 



(2.54) 



Due to the nature of the x{t) function (usually one chose a characteristic 
function), consequently in the expansion (12.501) we will take in account only 
the finite or singular part for t — > 

Our strategy to compute the action is to use the relation f)2.37p . therefore 
after explicitly expressed A and A\ we proceed to the calculus of the traces 
and in the end we will identify the leading part of the action comparing the 
result with the expansions fl2.50p - fl2.5ip . We can identify the operators A and 
A\ appearing in the (12.371) as follow: 



An 



(2.55) 



/ 



2i(l + ft 2 )L*(^)V 







(«) 







2i(l + n 2 )L*(Bv)V 



(n) 



with 



(2.56) 



Va,4, = (p^+il+n^iid^+A^), V B)4> = 4>* ( f>+(l+tf)(id tl B»+B tl *A> 1 ) 

(2.57) 
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VjP ) are define as V? } = (<9 M - itt 2 M.(x^)) and using the fl2T24l) it is 
useful to compute the commutators: 

[Vj?U*(/)] = ^(^/) (2-58) 
[HiMf)\ = -(l + ^ 2 )(^(A/) + 2L,(^)Vf ) (2.59) 

[^,VH = ^(Q-X^ (2-60) 

Referring to the f)2.37p and fl2.58p - fl2.60p for the computation of the action are 
required the following commutators 

[V\2i(l + £l 2 )L^A»)Vf ) ] = 

-(1 + Q 2 ) 2 {2iL i< {AA lx )Vf ) + ML*(d v A tl )Vf ) Vf ) ) 

[V 2 , [V 2 , 2z(l + Q^L^A^V^]] = 

8i(l + fi 2 )%(<9 p <9^)Vf V^Vj^ + • ■ ■ (2.62) 

From the previous discussion we are allowed to split the traces in two parts a 
matrix trace and the continuous one, now we focus ourself on the matrix trace 
contributions. Considering the matrix term e -* nE 4i2 ; com ing from e~ tQT,i , it 
contributes non trivially only for the vacuum tr(e~ tv ) and for the first order 
contribution tr(— tA e~ w ), in the other cases it contributes only by its leading 
part I32. The contribution from the vacuum is: 

tr (V^ 412 ) = 2tr (e"* 5 *) = 2tr ( f] e ^X-bA) 

On the base of C 14 defined by \si, ■ ■ ■ , s 4 ) = (b\y 
{0, 1} we have: 

tr ( e -tfffi4i a ) = 32 cos h 4 (fit) (2.64) 

The first order term related to the matrix trace vanishes, in fact it takes the 
form: 

t 2 fitr(S 4 r^)(Fr + FD (2.65) 

From f)2.2ip we deduce that the matrix trace tr(E 4 r^,) is proportional to 5^ u 
which contracted to the antisymmetric operator F^ u vanishes. The others 
matrix traces required are: 

tx{Y%Y%) = 4(1 + n 2 ) 2 (6, p 6 ua - 5M (2.66) 

tr^ + fir M+4 )r 9 • i(r„ + nr v+4 )r 9 ) = 16(1 + n 2 ) 6^ (2.67) 
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Turning our self on the calculus of the functional trace, to simplify the notation 
we introduce the functions : 

T(f) = Tr i2(K4) (L,(/)e-^) 

T n {f) = Tr i2(K4) (L,(/)V^---V^e-^) 

%u{f) = Trvm (^(/)Vg>V<2e-^) (2.68) 

Now substituting into (12.371) the (I2.55j) - (j2.56l) and using the matrices trace 
computations we obtain for the A field: 

Tr(e-^i) 

16cosh 4 (tfi)tr(e-^ 2 ) _ tT(lQV A ^) - t%(32i(l + n 2 )A») 

+t r (iqv a >4 , * v M + i6(i + n 2 )D^ *w^ + 8(i + n a )i^Ff 

+32z(l + ft 2 )^*<9 At V Aj(/ , 



t 2 / 

+ -% (32i(l + n 2 )A» * Vaj + 32z(l + ft 2 ) * A" 

\ t 2 

-64(1 + ft 2 ) 2 ^ * SU") + - V (-64(1 + ft 2 ) 2 ^ * A") 

-| % u ( - 64(1 + n 2 ) 2 v Aj4> * <9M* - 64z(i + ft 2 ) 2 ^ * ^y A ^ 

+64(1 + ft 2 ) 3 A p * (5 PI *AA V + 2d p dPA v ) 

-64(1 + ft 2 ) 2 (V^ * A" * A" + A" ★ V A ,4, *A v + A p -kA v * V A j) 
-128i(l + ft 2 ) 3 (A" * (<9 P A M ) * A" + A p * A"{dpA v 



t 3 

--% vp (128(1 + ft 2 ) 3 ^ * 9" A" - 128i(l + ft 2 ) 3 ^ * A M * A") 
+|V (512^ 2 (l + ft 2 )(e- 1 )^A^A p ) 

+^ V PCT ( - 256(1 + n 2 ) 4 A" * d v d p A a + 512i(l + ft 2 ) 3 ^ * A" * <9M C 
+256^(1 + ft 2 ) 4 ^ * («9^ p ) * A" + 256(1 + ft 2 ) 4 ^ * * A p * A' 7 



+5 field contribution + 0(Vt) (2.69) 

The contributions for the B fields are obtained operating the following substi- 
tutions: 

{A, B„, F M t F£, Vk,* Vb,«, ^ D^} (2.70) 
The next step is the computation of the T, T^...^ % the position space kernel 



Tuu can be ignored. 
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of e tHi (x, y) is 



- in i' X) y) = [ - — ] e -f (cothCn^l^-alHtanhC^la+yl 2 ) (2.71) 



27rsinh(2tfi) J 



which is essentially the four dimensional Mehler kernel (66] with £1 = 
and |x| 2 = x^x^. For the integral kernel of L*(/)(a:, y) we get after a direct 
substitution of (12.231) and a variable change: 

(Mf)9) = Jd'yU -0^f(x + ie • fc)e«<**->) y(y) 

A (pji / dS/(/)e 2 ^ fc ' - 1 ^ + ^ - 1 ^ + ^' - 1 ^)) <?(y) 

(2.72) 

From which the kernel of the operator L*(f) is easily identified as: 

L*(f)(x,y) = p- J d ^ z ) e i(--v,e-^+y))+Mz,e-H X -y)) ( 2 .73) 

Using f)2.7ip and (12.731) we compute the trace using a change of variables 
u = x — y, v = x + y and performing a Gaussian integration: 

Tr L2(R4) (L*U)e- tHi ) 

d 4 xd'ye- tH \y,x)LM)(x,y) 

& 

(27rfl) 4 27r 2 (l + Q 2 ) 2 sinh 2 (2tfi) 

--(i^+J&y)^ - 1 ^- 2 ^ 



d zf(z) d u d v 



xe 4 



cosh 2 (2tfi) 

cosh 2 (2tfi) J V ' J 

= =— ^- / d 4 zf(z)e- n ^ llz] (2.74) 

2tt 2 (1 + fi 2 ) 2 sinh 2 (2tfi) J V ' V ' 

In the end, we have: 

T(/) = =— ^- / d A zf(z)e- Q ^ llzl (2.75) 

V ; 2vr 2 (l + ^ 2 ) 2 sinh 2 (2tfi) J KJ 17 
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Using the same change of variables and the kernels fl2.7ip - fj2.73D we can com- 
pute %: 



VJ) 



1 f d*xd 4 yL^f)(x,y) ( A - 2i(&- 1 )^y v ) e- tH *(y,x) 



l + tt 2 J * * yJjy ,t)J \dy» 

1 



(1 + ft 2 )(27rfl) 4 47r 2 sinh 2 (2tfi) 



d z d u d v f(z) 



xe 4 



tanh(fit) coth(f2t) 



\i{u,Q- 1 (v-2z)) 



x u^— cosh(fit) — v^— tanh(f^) — zO (6 



2 ■ v "" v ■" v ~ JtIU 

1 f ,a „, , o tanh(flt) I 1 2 



(1 + fi 2 ) 3 4vr 2 sinh 2 (2^) 

with 



| dS/(^)e- fii w i|z|2 S^(e,7 7 ,t) | f= 

//\ tanh(fit) I |2 
dV(^ e (2.76) 



1 [ ,A . - _fitanh(flt)| ,2 



_ g ^tanh(flt) g q2/ q -i^ / 5 5 



(2.77) 



Z M = -fitanh^)^- — — z„ (2.79) 



2z'ft 

iTfi 2 

Following the same procedure %, u (f) turns out to be: 

%u{f) 



(i zf(z)e 



(1 + fi 2 ) 4 vr 2 sin 2 (fit) 

C7 2 /" otanh(Dt)i |2 

y (j^ + Z^ v ) e~ n ^^ lzl (2.80) 



(l + fi 2 )% 2 sin 2 (nt) 



with 



2tanh(mJ 



(2.82) 
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For the % vp (f) and % vpa {f) we obtain: 
Tfj,vp\f) 

d zf(z)e 1 1 



(1 + fi 2 ) 5 vr 2 sin 2 (fit) 

x qvd, + y^io* + y Vf &j v, t) \t= v=0 

= -—3 — / drzf(z)e i+« 2 1 1 

(l + fi 2 )% 2 sin 2 (fit) 7 

W(/) (2.83) 

~ (l + fi 2 ) 6 7r 2 sin 2 (^) 

IMjNvp ~\~ NfiyZpZfj ~\- J\fppZyZ a -\- N^fyZyZp ~\~ ZpZyJ\fp(j + ZpZpJ\f Va 

+Z^Z a Mu P + ZpZ u Z p Z^ (2.84) 

Finally, we have all the ingredients required to compute the leading part of the 
action (I2.45P replacing all the traces into the (I2.37p . Using the trace property 
of the star product and the identities 

*W = \{~^ /} + \[^ f] = /} + id J (2-85) 

{zp,f*g}* = {zp,f}*9-if*dpg (2.86) 

we get after some manipulations: 

Tr(e"^ 2 ) = 2 cosh 4 (fit) 



vr 2 (i + fi 2 ) 2 J \t\y Y i + fi 2 

/- 4fi 2 ~ 

+ ($*<!> + T^{xS*xZ-xl*xS] 

2V i + fi 2 B y 2 v i + ^ 2 M 
i ^(i-fi 2 ) 2 (i-fi 2 )^ /77lA , _ , _ , _ 



+ 2V 2 6(l + fi 2 ) 2 J ( F ^*^ + ^**b 

+ 1±^! (/^ *D^ + DpJ* Dy) \ + 0(y/i) (2.87) 
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where 

X^^ + Ap, X? = X° = ^ (2.88) 



Using the Laurent expansion of coth 4 (t') = t'~ 4 + ft'" 2 + § + 0(t' 2 ) 
1 ' Sit' 3SP 45 



W / ^{~((**# + I ^(^**J-^**ff)) 



TT^l + fi 2 ) 1 

+ 2 V 2 6(l + 2 )V l ^ A + ^ 5 J 

^^S^^p^Di^+D^-kD^) \(z) + 0(Vi) (2.89) 



Comparing the previous expression to the expansion (I2.5ip and putting xo 



x(0) we are finally able to write the spectral action (12.451) as: 

? 4 X- 4 , 20 2 X-2 , 52 Xo 



8fJ 4 30 2 45 

Xo 



z< 



4ft 2 



1 + ^ 



I + 02 "I 

+—— (Drf * DH + * D"^) (z) + 0( X i) (2.90) 
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Or after some rearrangements: 
S[Va) ~ W + + 




27T 2 (l + fi 2 ) 2 J I V I W I J- (")- V- / ■"" 1 /"' 



, na ^*JCj-^J +2(l + n 2 )^0* J D^^(:)-hC)( Xl ) (2.91) 

We notice that Higgs mechanism introduces an extension of the standard Higgs 
potential in the commutative case, in fact the Higgs scalar field <fi and the X%, 
Xg fields are present together in the potential. In this way the gauge field 
takes part in the definition of the vacuum. Another important property of the 
action, considering the X A , Xg as independent, is the invariance under the 
translations: 

4>{x) ->■ 4>{x+a), X^(x) -> X^x+a), X%(x) X%(x+a), Xg(x) -> X M (x+a) 

(2.92) 

which in others </> 4 -renormalizable theory is broken. Beside, the action is in- 
variant under U(l) x £7(1) transformations: 

<fi^UA*(fi*UB, X — > Ua * X^ A * Ua, Xg — )■ Mg * Xg *Ub (2.93) 

In field theory the ground state can be defined through the minimum of the 
action, the relevant part of the (12.9ip for the minimization is: 




2\4 



6(1 + n 2 ) 2 

2 

A- 1 \ 

A 



)K* f a+f>K v ) 



Xo 

2 



vB . i?H X—l 



1 + tt 2 



X»*Xl - ^) + 2(1 + &)Drf*Dv<f\{z) + 0(^2)94) 



Where we have omitted the constant part and we have rescaled the coefficient 
in front of the integral. Considering the fields X%, X^ as fields variables 
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instead A 11 , we can state that each terms of the action is semi-positive 
defined, so in order to find the minimum it is sufficient to minimize them 
separately. There are the two possible minima for the field strength part and 
for the covariant derivative part: the trivial solution with and X%, X^ equal 
to the null fields and the solution with 0, X A , Xg proportional to the identity. 
In each cases both the field strength part and the covariant derivative part 
disappear. For the potential parts we have: 



V A = V B =[ — \ for = X» = X% = 

Xo 



4Q y 
1 + " Xo 



rr ( 2 | A 4 ^ 2 X-l 

V B = (au + 4— — —a B 



1 + ^ 2 " Xo J 

Referring to the second case and minimizing the potentials, the minimum 
seems to be for 



W cosol, X» M = X» m = l^xP^smaV (2.95) 



™ Xo ' AU m 2 V xo V 4^ 

However, the previous position is not allowed because the identity does not 
belong to the algebra under consideration. In general the non-triviality of the 
vacuum makes very difficult to explicit the vacuum configuration of the system 
in [7T] A. de Goursac, J.C Wallet, and R. Wulkenhaar, using the matrix base 
formalism, have found an expressions from vacuum solutions deriving them 
from the relevant solutions equations of motion. Although, the complexity of 
the vacuum configuration makes the perturbative approach very complicated, 
in order to conduct some investigation in the next chapters will be consider 
a non-perturbative approach using a discretized matrix model of the action 
(I2.94p obtained using a particular Moyal base in which the fields become ma- 
trices, the star product becomes the matrix multiplication and the integral 
turns in a matrix trace. In this setting the action reduces to 




(1 - Q 2 ) 4 
6(1 + Q 2 ) 2 



yF**F? + F**Fj?) 



A , vfj, X-l 



2 



+ { ^ + —X^X A -- 




+2(i + n^*^ («) + °(xi) ( 2 - 96 ) 
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The omitted factor for the finite matrix model of size N becomes constant so 
can be ignored. The minimum is obtained like before and formally is (I2.95P in 
this case the identity, of course, belongs to the matrix space. It is interesting to 
notice that the vacuum of the finite model, due to the Higgs field, is no longer 
invariant under the transformations (12. 93p . but is invariant under a subgroup 



Having discretized the model will be performed a Monte Carlo simulation 
studying some statistical quantity such the energy density, specific heat, vary- 
ing the parameters Q, a and gathering some informations on the various 
contributions of the fields to the action. The simulations are quite cumbersome 
due to the complexity of the action and to the number of independent matrix 
to handle. Anyway, using particular algorithm and some simplifications we 
are able to get an acceptable balance between the computation precision and 
the computation time. 



of U(N) x U(N): 




(2.97) 
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Introduction to numerical 
analysis 



In this chapter will be explained the basics of the Monte Carlo simulation, 
focusing ourself on the application of such simulation on the field theory. In 
the first part will be hinted the path integral formulation which is essential 
to connect the field theories to statistical systems in which the Monte Carlo 
methods are born. Then will be introduced the Metropolis algorithm used to 
produce a Markov chain. In order to resume all the previous concepts and to 
show an example of phase analysis will be presented an example of numerical 
simulation on the well know Ising model. 

3.1 Path integrals and functional integrals 

The path integral formalism in quantum mechanic was introduced by R. P. 
Feynman as a generalization [621 |6T] and provides a powerful tool to study 
quantum field theories. The path integral is based on the superposition law, 
in fact one has to consider a superposition of all possible paths in order to 
compute the transition amplitude from an initial state at to to a final at time 
t\. Considering the initial and final states written as \ip(x )) at t and \ip(xi)) 
at time ti, the transition amplitude is determined taking the matrix element 
of the temporal evolution: 



is the time evolution and H is the Hamiltonian of the system, which we assume 
to be time independent. The time interval (t ,tx) can be divided into N 



(V>(zi) | U{t u t ) | V(aro)) 



(3.1) 



where 




(3.2) 
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sub intervals of size e = r i+1 — r i7 with t = r and ti = t n , in this way is 
possible to write the time evolution operator as: 

-iH!izio) = p _^ (r ^_ Tjv _ 1+Tjv _ 1 _..._ T2+r2 _ Tl+Tl _ To) = f p -^y (33) 

2 

In the case H has the general form H = Hq + V, with Hq = we can write 
using Trotter formula 

for e — )■ and supposing that ifo and are semi-bounded. In order to factorize 
(13. ip is useful to consider a set of complete states 

(xi)){t/)(xi)\dxi = 1, x 1 = x(Ti) (3.5) 

and inserting it between each term e - tH o e / h e - tVe / h . As result the (13.11) can be 
rewritten as N terms product 



{tp(xi) | U(t 1 ,t ) | ip(x )) ~ / dxAr_i / dxAT_ 2 -- - / dx 



i 



x(<f)(x N ) I e ~«be/« e -iV e /R | 1 /,(5 JV _ 1 )) 

X(0(x o ) I e - iHoe ^e- iVe / n | (3.6) 

Considering a conservative l^(x) and H dependent only on p, it is possible to 
demonstrate that EH [63] : 



(3.7) 

In the previous statement the argument in the exponential can be written as 

^ X ' - H(p i+1 , ^(x i+1 + Xi)j (3.8) 

It easy to recognize the action in discrete interval (r^, Tj+i) like the Lagrangian 
in the interval times the duration of the interval. Using ( 13. 7p and ( 13. 6 p taking 
JV^oowe finally obtain: 

(iP( Xl ) | t ) I V(*o)> = f (M 11 e- 5 W/ fi Dx(t) (3.9) 

Jx(to)=X 
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where 

S[x] = / L(x,x)dt, (3.10) 

J to 

S[x] and L are the action and the Lagrangian of the system, Dx(t) stands 
for the integral over all paths x(t) starting from x(£q) = xo to x{t\) = x\ and 
Dx(t) is the functional measure. The generalization (13. 9p to quantum fields 
is straightforward. For completeness, before introducing its expression, it is 
important to remind that in the path-integral methods, when Osterwalder- 
Schrader axioms hold [BHE5], it is common to have the action with imaginary 
time in order to simplify the calculations. Beside, there is a strong advantage of 
considering an imaginary time, in fact it allows to establish a relation between 
statistical mechanic and quantum field theory. 

The euclidean frame approach is based on the analytic continuation, which 
is a technique used in the domain of definition of a given analytic function in 
order to extend a real valued function into the complex plane. Our purpose is 
to start from real time t to the imaginary time t, the so called the Euclidean 
time [63] . In this way we force the imaginary time and spatial coordinates to 
have the same signature. A simple application is the D'Alembertian operator 
which in real time is given by: 

□ _ na_ _ _ _ an) 

dt 2 dx\ dx\ dx\ 
Applying the previous procedure we obtain: 

A - - — (312) 

dxl dx\ dx\ dx\ 

with xq = t = it. In quantum mechanic, when we consider the path integrals 
over all possible particle paths between two points, we treat space and time 
in different ways. Otherwise, in field theory, in order to treat them in the 
same manner, we introduce a function of 4-dimensional space-time manifold 
x = (x,t), which we call field ip(x). In the case ip is a neutral scalar field, the 
field values are real, ip(x) G R. We define a configuration like a fixed field value 
in each space-time point x, the configuration plays the role of particle path 
in quantum mechanic. The functional integral is over all field configurations 
J Dip, the Lagrangian L becomes a Lagrangian density L(^(x),d /I i/}(x)) at 
each point x and the action is obtained by an integral over the space-time 
volume S = J Ld x. The generalization of (13. 9p for quantum fields is given by 
(h = l): 

(AlTihMA) = [ e- s ^(x,t) (3.13) 



Beside the path integral there is also a second interpretation of the resulting 
functional integral, as a functional integral in statistical field theory. In this 
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frame the Euclidean action is seen as the energy functional of an analog sta- 
tistical mechanical system with k^T = 1 where ks is Boltzmann's constant. 
A statistical system in equilibrium is naturally described by ( I3.13p . 

3.1.1 Expectation values 

The expectation values of an observable O it is defined as follows: 



is called partition function. The analytically integration in (I3.14p . in general, is 
impossible, indeed it involves all the possible configurations in the functional 
space. But we can try, using numerical methods, to estimate the value in 
f)3.14p . One of the most popular approach is the Monte Carlo method, this 
method is founded of the assumption that we can estimate the expectation 
values using representative sets of random configurations. There are many 
way to generate a representative set, usually are used random moves from 
a starting configuration in order to explore the configuration space. In this 
thesis, to estimate the expectation values of the observables defined in the next 
chapter, will be used a particular algorithm for the configurations generation 
called the Metropolis algorithm |75j . 



Many physical and mathematical systems are often treated using Monte Carlo 
methods thanks to their computational algorithms which use a random (or 
pseudo-random) generated sampling, these methods are most suited to calcu- 
lation by a computer and tend to be used when it is hard or impossible to 
compute an exact result with a deterministic algorithm. These methods are 
widely used in mathematics; one of the classic application is the evaluation 
of definite integrals, particularly multidimensional integrals with complicated 
boundary conditions. Beside, Monte Carlo methods have a great importance 
in statistical mechanic for example in the statistical systems with a large num- 
ber of degrees of freedom, such as strongly coupled solids, disordered materials 
and fluids. A precise definition of Monte Carlo methods is hard to give, in fact 
there is no single Monte Carlo method, instead the term describes a large and 
heavily used class of approaches and algorithms. However, these approaches 




(3.14) 



where 



Z 




3.2 Monte Carlo methods 
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tend to follow a particular procedure. The first step is to define a domain 
of possible inputs, generate inputs randomly from the domain using a certain 
specified probability distribution, then perform a deterministic computation 
using the inputs and at the end join the results of the individual computations 
into a final result. For the present thesis one possible definition Monte Carlo 
could be: a method of approximating an expectation value by the sampled 
mean of a function of random variables, invoking the laws of large numbers. 

We can formalize the previous ideas [75] considering a random variable 
x having probability function or probability density function px(%) which is 
greater than zero on a set X of values. Then the expected value of a function 
/ of the continues variable x is 

E(f(z)) = [ f(x) Px (x)dx (3.15) 

For the discrete case : 

E(f(X)) = ^f(x)p x (x) (3.16) 

A Monte Carlo estimate of (13.15P or (I3.16P can be defined taking N samples 
of X, (xi, ■ ■ • ,Xn) and computing the mean of f(x) on the sample: 

1 N 
1=1 

Jn{x) it is called Monte Carlo estimator of E(f(X)). This estimations can be 
applied either when the generated variables are mutually independent or when 
they are correlated one to another (for example if they are generated by an 
ergodic Markov chain). For simplicity we will consider the random variables 
independent, but all can be extended to samples obtained from a Markov 
chains via the weak law of large numbers. If E(f(X)) exists, then the weak 
law of large numbers tells us that for any arbitrarily small e: 

lim P(\f N (x) - E(f(x))\ > e) = (3.18) 

JV— >oo 

In other words, for N large the probability that /jv(^) deviates much from 
E(f(x)) becomes small so we are justified to use (I3.17P as estimator of E(f(x)). 
It is interesting to note that expectation value of the (I3.17P is: 

/ 1 N \ 1 N 
E(f N (x)) =E[-Y, ) = N £ WW = E M*)) (3-19) 

\ i=l / »=1 

In this case we say that /jv(^i) is unbiased for E(f(x)). The previous method 
turns to be very useful when it is applied in all the situations in which the 
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quantities of interest are formulated as expectations value, for example: prob- 
abilities, integrals, and summations. 

Let Y be a random variable, the probability that Y takes on some value in a 
set X can be expressed as an expectation using the function: 

P(Y e X) = E(C (X) (Y)) (3.20) 

where C X (Y) is the characteristic function of X that takes the value 1 when 
Y E X otherwise is 0. 

For integrals we consider a problem which is completely deterministic the 
integration of a function f(x) from a to b. This integral can be expressed as 
an expectation respect to a uniformly distributed continuous random variable 
u between a and b, with density function P u {u) = 1/(6 — a). Rewriting the 
integral we obtain 

(b-a) I ^-dx = {b-a) [ q(x) Pw u(x)dx = (b - a)E(q(u)) (3.21) 
Jb °-a J b 

A discrete sum is just the discrete version of the previous example; the sum 
of a function q(x) over some numerable values of x in a set X. Using a 
random variable w which takes values in X all with constant probability P 
with Ylwex P — ^ the sum can be seen as the expectation: 

E = JS E 9(*) p = 7^(9(«0) (3-22) 

The immediate consequence is that all probabilities, integrals, and summa- 
tions can be approximated by the Monte Carlo method. However, it is very 
important to point out that there is no a priori reason to use uniform distri- 
butions. As we have seen many quantity of interest can be formulated as an 
expectation approximate by a Monte Carlo estimator, but it is not always so 
easy to actually have a Monte Carlo estimator that can provide a sufficient 
good estimation in a reasonable amount of computer time. For the same prob- 
lems various number of Monte Carlo estimators can be constructed (essentially 
varying the probability distribution), of course some Monte Carlo estimators 
are more efficient than others. In order to find the best estimator we need to 
compute the variance, we look for the Monte Carlo estimator with smallest 
variance taking the amount of computational effort fixed. We have to com- 
pute the variance of the Monte Carlo estimator of E(f(x)) with the random 
variable /at(x). The standard formulas for a random variables are 

Var(MX)) = VorUt /«)) = 

= El/W-TOlf^W (3-23) 

xex 
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if x is discrete, and 

= / [/(x)-B(/(A-))] 2 P v (x) (3.24) 

JxeX 

if x is continuous. There are numerous sophisticated algorithm used to obtain a 
better approximation reducing the variance, many of them uses a non-uniform 
probability distributions, in the next section we introduce a such kind of algo- 
rithm called Metropolis algorithm which is particularly useful for the systems 
with many degrees of freedom. 

3.2.1 The Metropolis algorithm 

[73] Studying systems with a great number of degrees of freedom it becomes 
clear that the number of possible configurations becomes exceedingly large very 
quickly. Even for extremely simple binary models, in which each particle may 
exist in one of two possible states, the number of configurations of the system 
grows extremely rapidly with N, the number of particles. A very instructive 
and simple example is the Ising model of magnetic spins, in which each spin- 
1/2 particle can be spin up or spin down respect to a fixed axis. In this way 
each of the N particles has only two possible states, so the total number of 
discernible configurations is 

In easy to notice that using a model of a square lattice of such spins taking a 
quite small scale of a macroscopic sample of matter with 32 x 32 = 2 10 spins, 
there are 

2 2l ° w 10 308 

distinct configurations of this system. But remind the basic postulate of sta- 
tistical mechanic, we know that an isolated system in equilibrium is equally 
likely to be found in any one of these configurations. To compute the time- 
averaged properties we compute an ensemble average over all the accessible 
configurations. If system is not isolated from its external environment, but 
may exchange energy with it, the exchange is ruled by a temperature T, which 
quantifies how the energy is shared with its surroundings. The Boltzmann fac- 
tor e~ E l T is proportional to the probability that the system will be found in a 
particular configuration at energy E when the temperature of the environment 
isT(K b = l) 

P oc e~ E/T (3.25) 
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Higher T means that the environment tends to give energy to the system 
otherwise in the small the temperature case the environment tends to gather 
energy bringing the system in a lower energy state. In view the identification of 
(I3.13P as description of a statistical system in equilibrium (fc^T = 1), explained 
in the previous section, we can repeat the same treatment in the case we have 

P oc e~ s (3.26) 

In the study of a such model, unless we can treat it analytically, it is very 
hard to try investigating each configuration, or to average over all of them. 
Usually when we handle an interacting system, we can almost never compute 
the sums analytically, so we have no choice but to seek an approximation 
either analytically or numerically. One choice is to use the a Monte Carlo 
approach with the Metropolis algorithm in order to limit the huge number of 
possible configurations. The Metropolis algorithm is based on the notion of 
detailed balance that describes equilibrium for systems whose configurations 
have probability proportional to the Boltzmann factor. Essentially we sample 
the space of possible configurations in a way that agrees with (13.26P . Beside we 
complete the exploration by taking in account the possible transitions between 
close configurations. Consider two configurations A and B, each of which 
occurs with probability proportional to (I3.26P and computing the ratio of the 
probability of the configurations we obtain 

P (^) _ e ~ SA - P -(S A -S B ) (o 27) 

The notable thing about considering the ratio is that it tuns a relative prob- 
abilities involving an unknown proportionality constant (the inverse of the 
partition function) into a pure number. Metropolis et al. [75] have proposed a 
practical algorithm to achieve the relative probability of (I3.27p . In a simulation 
the procedure is the following: 

1. Beginning from a configuration A, with known Sa, make a change (typ- 
ically small) in the configuration to obtain a new configuration B. 

2. Compute Sb (usually do not differs so much from Sa)- 

3. If Sb < Sa, accept the new configuration, since it has lower energy 
according to (13. 26 p . 

4. Otherwise the number e~^ Ss ~ SA ^ is compared to a random number, ran e 
[0,1]; if e~( Sj3 ~ SA ^ > ran the configuration B is accepted, in the other 
case it is rejected. 

5. If the configuration B is accepted go to point 1 using B as new starting 
configuration, in case of rejection go to point 1 trying another configu- 
ration B. 
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Following these prescriptions, we will sample points in the space of all possible 
configurations with probability proportional to the exponential factor consis- 
tently with the theory of statistical mechanic equilibrium. We can approximate 
expectation value like (I3.17P by summing them along the path generated by 
Metropolis procedure. Some hard parts about implementing the Metropolis 
algorithm are the first step, how to decide the staring configuration and how 
to generate useful new configurations, the first issue will been solved using the 
thermalization graphs. 

3.2.2 Phase analysis 

[79] In order to show and test advanced analysis tools for stochastically gen- 
erated data it is convenient to work with very simple models, in fact a very 
few exact analytical results are available for comparison. On the other hand, 
one should always be prepared that a really complex system may add further 
complications which are not present in the simple test cases. Nevertheless, 
to introduce the phase transitions we go back to the well known Ising model 
whose partition function is defined by the Hamiltonian Hf. 



where (3 = J/ksT is the inverse temperature in natural units, the spins take 
place on the sites i of a n-dimensional cubic lattice of volume V = L n and 
the symbol (ij) indicates that the lattice sum is on all nearest pairs assuming 
periodic boundary conditions. In two dimensions and with no external field 
this model has been solved exactly even on finite lattices. Unfortunately for the 
three dimensional model is not so easy and an exact solution is not available. 
Instead, there is an huge amount of very precise Monte Carlo simulations. 
The standard observables are the internal energy density e = E/V, with E = 
d\nZ I /d(3(H I ), the specific heat 




(3.28) 




(3.29) 



the magnetization density 




(3.30) 



and the susceptibility 



(H> 2 ) 



(3.31) 



Beside, the spin-spin correlation function 



can be defined as 




(3.32) 
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Can be showed that at large distances, the spin-spin correlator G(x) decays 
exponentially G(x) ~ e - ' 2 ^, so we can define the so called spatial correlation 
length £ as 

£= lim -^0— (3.33) 
\x\^oo kiG(x) 

The theory of phase transitions is a large subject, here we shall confine ourself 
to those properties that are important for understanding the basic concepts 
for the data analysis. Most phase transitions in nature are of first order and 
are characterized by discontinuities in the order parameter Am of the mag- 
netization m in Fig. 1 or in the energy like latent heat Ae or both, at the 
transition point T . This is connected to the coexistence of two or more phases 
at To. In contrast to a second-order transition, the correlation length in the 
coexisting pure phases is finite. Therefore, the specific heat, the susceptibility 
and also the autocorrelation time do not diverge in the pure phases. However, 
can be present superimposed delta function like singularities associated with 
the jumps of e and m. The standard example exhibiting a first-order phase 
transition is the g-state Potts model defined by the Hamiltonian 

H P = -J2^ , ^el,--- ,q (3.34) 

(id) 

where 8ij is the Kronecker symbol. This model in two dimension is exactly 
solved and exhibits a first order transition at (3 t = log(l + y/q) for all q > 5 
and for q < 4 the transition is of second order, including the Ising case (q = 2) 
and the special percolation limit (q — 1). In three dimension there are many 
numerical evidences that suggest that for q > 3 the transition is of first order. 

The characteristic feature of second-order phase transitions is a divergent 
spatial correlation length £ at the transition point /3 C . This causes a scale 
invariance, the core of renormalization group treatments [77] and is the origin 
of the important concept of universality. At (3 C are expected fluctuations on 
all length scales, implying singularities in statistical functions such as the 
correlation length 

£ = + (3.35) 

where t = |1 — T/T c \ and the • • • indicate higher order corrections. In the 
previous expansion are defined the critical exponent v and the critical ampli- 
tudes £j) where ± symbols denotes the high and low temperature side of the 
transition. Singularities are present in the specific heat, magnetization and 
susceptibility too fig J3.ll [79J. They define, respectively, the critical exponents 
a, (3 and 7, can be proved that these coefficients are related each other through 
scaling relations and only two of them are independent. 

When the spins are updated with an Metropolis-Monte Carlo process, [79] 
the information on the updated state of the spins has to propagate over the 
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correlation volume before one obtains a new truly independent configuration. 
The number of update steps required is measured by the autocorrelation time 
r (its definition is given in appendix) which close to (3 C behaves according to 

t oc C oc r vz (3.36) 

Here we have introduced the independent dynamical critical exponent z but 
depends on the used update algorithm. For the Metropolis we expect that the 
updated information performs a random walk in configuration space, requiring 
on the average £ 2 steps to propagate over a distance proportional to £ therefore 
we expect a dynamical critical exponent of z ~ 2. Of course in systems 



Ftrsl -Srdffl TraisfijQfl Ssoand -Order Ttan&lra n 





Figure 3.1: The characteristic behavior of the magnetization, m, specific heat 
C and susceptibility, x densities at Grst and second order phase transitions. 

characterized by finite size, as in any numerical simulation, every statistical 
function like the correlation length cannot diverge and the divergences in all 
other quantities occurring in the infinite model are rounded and can be shifted. 
The example of the specific heat density for the two dimensional Ising model 
is illustrated in fig l3.2l The finiteness affects the scaling formulas and £ is then 
substituted by the linear size of the system L. If we write 

t oc C 1/v -> L- 1/u (3.37) 
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it easy to understood that the thermodynamic scaling laws x ^ > f° r finite 
models, has to be replaced by finite-size scaling (FFS) x °^ L 1 ^ . In particular, 
[77] we obtain for the autocorrelation time a FSS of the form 

r oc L z (3.38) 



«2 
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Figure 3.2: Finite size scaling behavior of the specific heat density of the 2D 
Ising model on L x L lattices close to the critical point j3 c . 

Even for first order phase transition in a finite systems the singularities can 
be split in narrow peaks with height proportional to the characteristic volume 
V and width proportional to V~ x . The system is now capable to exist in a 
mixed phase configurations and jump from one phase into the other, however 
the probability of the appearing of this mixed configurations is smoothed by 
the Boltzmann factor: exp(2crL™~ 1 ). The double-peak structure of the energy 
and magnetization densities that can occur for a first-order phase transition 
too and can be explained by the previous argumentations. For large system 
sizes it may require many Monte Carlo steps before the systems explore enough 
phase space. In fact, the autocorrelation time associated with this transitions 
is the inverse of the suppression factor: 

t oc exp(2aL"" 1 ) (3.39) 

Since in this case the autocorrelations grows exponentially with the system 
size, this behavior is usually called supercritical slowing down. 
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In this chapter will be set up all the ingredients required to implement a Monte 
Carlo simulation of the model treated in the chapter two. The first part will 
be devoted to the discretization scheme, introducing a Moyal base in which the 
model described by the action (I2.96P on Moyal plane will be discretized. In 
the second part, in order to define the measurements of the conducted Monte 
Carlo simulation, will be applied some arguments treated in the third chapter 
outlining the expectation values, statistical quantities and some order parame- 
ters. 

4.1 Discretization of the action 

The first step across the numerical analysis is to apply a discretization scheme. 
Various schemes can be used like lattice approximation, but the nature of star 
product due to its oscillator exponential, makes the lattice approach not suit- 
able without adaptations. We will use another discretization scheme in which 
our fields are approximated by finite matrices and the star product becomes 
the standard matrix multiplication. Our model is Euclidean, as remarked in 
chapter 3, the main advantage of working in this formalism is that we can 
establish a direct connection to statistical physic and the functional integral 
converges with a relatively modest statistic. The action to discretize is 



S(4>,X A ,X B ) 



I 




(i-» 2 ) 2 (1-r; 2 ) 4 

2 6(1 + ft 2 ) 2 




> r /MU * r B 




(4.1) 
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Using the identity 
following form: 



S(<P,Xa,X b ) 



d 4 z< 



i-n 2 



Xa/j, * we can recast the action in the 



i-n 2 



6{i + n 2 y 



Xbh, Xbv 



An 2 



^B'^B 



X^-kX 



Xa/j,, Xa 
40 2 



v-M yv 



X A * Xau 



i + n 2 B 



By, 



X-i 
Xo 



+2(1 + n 2 ) (V * x Bil - x Afi * 4>) (« 



i + n 2 A 



X-i 
Xo 



+o(xi) 



(4.2) 



As a first approach to the numerical simulation and forced by limited compu- 
tation resource, we will consider the Monte Carlo simulation of the previous 
action around its minimum and the simulation will take \J^~ as a positive 

parameter. In this setting the behavior of the simulations will be identical for 
the negative case and avoiding - ^- to be negative we have not any problems 
about the thermalization. In order to define the previous action around the 
minimum we translate the fields 0, Xa/j,, X Bf i using the following translated 
fields: 



cosal 



X 



An 



X 



tin 



Substituting the previous fields into 
mum in zero: 




2Q 2 



2Q 2 



(4.3) 
(4.4) 

(4.5) 



we get a positive action with mini- 



s(^,y a ,y b ) = yV*|i?([K^,y^^ 

+ (V> * $ + fi cos a(ip + 4>) + CY^ ★ Y Ail + I^YY Afl sin a) 2 
+ (ip*^ + H cos a(if> + + CYg * Y Bfl + fiYY Bfl sin a) 2 

+2(1 + Q 2 ) {{Y Blx - Y Alx )n cos a + V * Y Bll - Y A[X * i/>) 

*((y^-y^cosa + ^*y^-y^*^) J(z) + e>( Xl ) (4.6) 

64 



Chapter 4. Preparing for numerics 



Where for simplicity we put: 

Recalling the first chapter we consider a Moyal base in which our fields became 
infinite matrices and the integral becomes a double trace, we will discretize 
the action introducing a cut off and then we will show how to simplify the 
discretization in order to conduct a feasible numerical simulations. 



4.1.1 Discretization by Moyal base 

The following treatment is mainly taken from [55j [29] as introduction to the 

Moyal base which will be used later. We can define on the algebra M a natural 
basis of eigenfunctions f mn of the harmonic oscillator, where m, n G N. If we 
define the Hamiltonian 

H = H 1 + H 2 + --- + H N , (4.8) 

with 

H r .= ^ + xl N )iorl = l,--- ,N (4.9) 
then the f mn diagonalize these harmonic oscillator Hamiltonians: 

H -k f mn = 6{m + -)f mn , fmn * Hi = 9(n + -)f mn (4.10) 

In order to derive f mn we note that the Gaussian function 

fo(x) = 2e~^ +x ^ (4.11) 

is an idempotent, 



(/o*/o)(a:) 



We can define the graded creation and annihilation operators 

1 (9 1 

a = -j=(x 1 -ix 2 ), — = + ify) (4.13) 

Computing the left and right star multiplication of this graded operator with 
any / G R| we obtain 

9 df 9 df 

(a */) = a(x)f(x) + 2~fa( x ^ (f * a ) = a ( x )f( x ) ~ g -q= ( x ) 



9 df 9 df 

(a ★ /) = a(x)f(x) ~2~fa( x ^ (/* a ) = H x )f( x ) + 2 



(4.14) 
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The multiple application of the star multiplication for a and a to /o leads to 
a* n * f = 2 m a m / , fo * = 2 n a"/ (4.15) 

and 



a * a* m -k /o 



m^(a*( m " 1 ) * /o) for m > 1 
for m = 



^(/o^a*^- 1 )) form>l 
/ *a*"*a = <J (4.16) 

for m = 

with a* n = a* a* ■ ■ ■ -k a and similarly for a* n = a * a * • • • * a, where in the 
left side appear n factors. Finally, can be proved by induction the following 
relation for f mn : 

f mn = -. 1 a* m * f * a* n (4.17) 

min(m,n) 

■\m-\-n—2knk — m—kn—k 



— r (-i) fe f?V?V'2 r 



c nk — m—k „n—k f 

v x a a Jo 



\/n\m\i 1 k=Q 

From the equations f)4.12p and (I4.16P we obtain the ^-multiplication rule: 

(fmn * /fcj) (a?) = 5 nfe / m J (x) (4. 18) 

The previous multiplication rule associates the ★-product between f mi with 
the ordinary matrix product: 

oo oo 

(l(x) ^ ^ Q"mnfmn {•£ ) j b{x) ^ ^ b mn f mn (x^j 
m,n=0 m,n=0 
oo oo 

(d*&)(x)= 2j (^)mn/mn(l), Mm n = ^]flmAn (4.19) 
m,n=0 fc=0 

In this base we can write any elements of 1R@ but we have to require the rapid 
decay [29] of the sequences of coefficients {a mn } appearing in the (14.191) : 

oo oo 

a mnfmn(x) E R|, if, ^ (( 2m + + ^Knl^ < OO for all k. 

m,n=0 m,n=0 

(4.20) 
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Using again (14.121) . the trace property of the integral and (14.161) we find this 
useful property: 

d 2 xf mn (x) = - 1 / d 2 x{a* n * f Q * f * a* n (z)) = 5 m „ / d 2 xf 

^n\m\6™ + J J 

= 27c65 mn (4.21) 
The eigenfunctions f nm can be expressed with the help of Laguerre functions 



n—m 

f mn (P,v) = 2(-l) m ^e^- m ) f ^p) e^L— Qp 2 ) (4.22) 
Our fields can be expanded in this base as: 



X"(x)= V" X^ ini f mini {x ,xi)f m2n2 (x 2 ,x 3 ) (4.23) 

' ■ m 2 n 2 

and 

V'(aj) = ^^i^/mxm^O^O/mana^,^) (4.24) 

Using this base we can forget the Moyal product in this way the model becomes 
to 9-matrix model. In view of (14.191) a ^-product between two fields using 
(I4.18P can be written as 

= } qfm ini ® klh f (x ,Xl) */ Wl (x ,Xl) 

^ ' " l 2™2 k 2 l 2 

mi,ni,fci,ii£N 
X fm 2 n 2 (x 2 ,X 3 ) -k f k2h (x 2 ,X 3 ) 

= y ^ m i"i^k 1 i 1 5 nikl S n2k2 f mi i 1 (x ,x 1 )f m2 i 2 (x 2 ,x 3 ) 

'™2 n 2 k 2 l 2 

mi,ni,fci,iiSN 

V" ^§m 1 i 1 f mih (xQ,Xi)f m2l2 {x 2) x 3 ) (4.25) 



rrti,ZieN 



where 



rr, n U rn 9 no „„;„ 



m 2 Z 2 m 2"2 n 2 i 2 

ni,n 2 eN 



So the star product became a "double" matrix multiplication, the action, 
the equations of field and all treatments can be conducted on the infinite 
matrices instead directly on the continues fields. Beside, for finite matrices, 
the N 2 -indexed double sequences can be written as tensor products of ordinary 
matrices, 

K 



Xm 1 n 1 =S^X i 0X i . (4.27) 

mono. /-^ m l"-l ^ m 2 n 2 ■ K^-^'J 
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Since the matrix product and trace also factor into these independent compo- 
nents, the action factors into S = "£f =1 S(ip u , Y%, Y^)S(ip 2i , Y 2 \ Yf). Then, 
regarding all iI> 1i ,Ya, Yg ,ip 21 ,Y^ ,Yg as random variables over which to in- 
tegrate in the partition function, the partition function factors, too: 



J Y\\ Yh\^\ Y 2 \ Fj 1 ) ■ ■ ■ V(^ K , Yl K , Yl 

V^ l \Y l A \Y l B l ^ 2 \Y 2 \Y 2 ') e -sW\Y%,Y%).sw>,Y X >,Y S n ) ( L2(S) 



We may therefore restrict ourselves to K = 1 and we can split the model in 
2+2 dimensions, in this way the calculus will be performed just on standard 
matrices. This decomposition makes the numerical treatment of 2 dimensional 
case simpler respect the complete one, indeed will be examined both cases. 
Using this base now our problem is reduced to a infinite matrix problem but 
is not enough to be handled numerically, we have to perform a truncation 
in order to obtain finite matrices, this truncation will consist in a maximum 
m,n < N in the expansion (I4.24p - fl4.23p . It is easy to verify that this kind of 
approximation corresponds in a cut in energy, in fact recalling the definition 
f)4.10p of j mn we have: 

{H, f mn }* = 6(m + n + ~)/ mn (^ p) (4.29) 

Beside, the functions f mn with m,n < N induce a cut-off in position space 
and momentum space, indeed the function 

L^(z)z a/2 e- z/2 (4.30) 

with z = |p 2 , is rapidly decreasing beyond the last maximum {z^) max . Nu- 
merically can be found (z^) max < 2a + 4m and thus the radial cut-off 

p max ~ V29 , for m,n < N (4.31) 
In addition, for p 1 = psinip, p 2 = pcoatp we obtain [55] using f )4.22p , 

/>00 />27T 

/m»(p,V0 = / Pdp j PpSiu{ *-^f m n{pM 

Jo Jo 

/ i — \ n—m 



/ , — \ n — rn 

i(*+*)(n-m) U e -p\ L ^(t p ^ e ~iP 2 (4.32) 



n — m 

, / m' 
2n8 



In the end we obtain 



Pmax ~ \l —pr for m,n < N (4.33) 
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Summarizing, to operate the discretization we have the following correspon- 
dences: 

<p(x) e Kq -> 4> e M N (4.34) 
#£»e -> y/6M ff (4.35) 
^WeRe ?eMi (4.36) 

Tr(a) (4 - 37) 

After truncating the representative matrices is convenient to operate another 
substitution I7T1: 



Zo 




+ iYf, 


Zo 


- Y A 


— iY A 






+ iY*, 


Zi 


- Y B 


— %Y B 


z 2 


- Y A 


+ iY A , 


z 2 


- Y A 

— I 2 


— iY A 


z 3 


- Y B 


+ i% B , 


z, 


- Y B 





2 



2 



(4.38) 

Implementing the previous replacements for the two dimensional model, thus 
using only Z , Zi, in eq. (14.61) we arrive at 

S 2 = Tr (£p + C Vo + C Vl + C Do C Do + C Dl C Dl ) (4.39) 

with 

C 2F = ^([Z ,Z ] 2 + [Z 1 ,Z 1 ] 2 ) 

C 2Vo = ^ + ^coaa(^ + ^) + ^{Z ,Z } + ^^^((-l + i)Zo + (l + i)Z )^ 

C 2Vl = ^ + ficosa(i> + 4>) + ^{Z 1 ,Z 1 } + ^^{{-l + i)Z 1 + (l + i)Z 1 ) 

C 2Do = y/2(L + ( M cosa(Zi + Z x - Z Q - Z Q ) + ${Z X + Z x ) - (Z + Z )i/>) 
£ 2Dl = ^2(1 + ft 2 ) (/icosa(Z 1 -Z 1 -Z Q + Z Q ) + ^{Z x - Z x ) - (Z - Z )i/>) 

Where for simplicity we have omitted the hat on the matrices and the bars 
stand for the hermitian conjugate. With this discretized action (14. 39p . the 
scalar if> and the gauge fields Zq, Z\ become 3-complex matrix instead to have 
one complex matrix for the scalar and four real matrices for the Y^, Yb^ fields. 
In this setting the algorithm for Monte Carlo simulation becomes simpler es- 
pecially with the update algorithm used here. In the four dimensional model 
it is convenient to use the (I4.38[) too, in this case there are some additional 
terms. In particular appear the crossed commutators between the Zi in the 
Y-M part: 

Si = Tr [Cp + £y Q + Cvx + £d £-D + £-D 1 £-D 1 + £-D 2 £-D 2 + £d 3 £d 3 ) (4.40) 
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With 





= ^[z Q ,z a ] 2 +[z 1 ,z 1 ] 2 + ^[z a + z ,z 2 -z 2 _ 


2 -[ 


-Zo + -^o, Z 2 + 




+ [Zq — Zq, Z 2 + Z 2 ] — [Zq — Zq, Z2 — Z 2 ] — [Zi - 


f z u 


z 3 + z 3 ] 2 




+ [Zi + Z U Z 3 - Z 3 ] 2 + [Zx - Z 1: Z 3 + Z 3 ] 2 ~ [Zx - 


-Zx,Z 3 -Z 3 ] 2 ^ 




= (^ + /xcosa(^ + $ + - {{Z ,Z } + {Z 2 ,Z 2 }) 
/isina 1 + t )( Zo + Z 2 ) + (1 + i)(Z + Z 2 ))) 2 

2 V L/ 








= $i/> + M cosa(^ + i/0 + i ({Zi,Zi} + {Z 3 , Z 3 }) 
+ ^-((-1 + 0(^x + 30 + (1 + + Z 3 ))f 






£iD 


= a/2(1 + ft 2 ) (^cosa(Zi + Zi - Z Q - Zq) + ip{Z x 4 


-^1) 


- (Z Q + Z )if>) 




= ^2(1 + Q 2 ) (^cosa(Zi - Z 1 - Z + Zq) + ip(Zi - 




- (Z Q - Z )fl>) 


£-4D 2 


= y/2(l + (vcosa(Z 3 + Z 3 -Z 2 - Z 2 ) + i>(Z 3 4 


-^3) 


- {Z 2 + Z 2 )if>) 



U Ds = y/2(l + W) (MCOsa(Z 3 -Z x -Z 2 + Z 2 ) + ^(Z 3 - Z 3 ) - (Z 2 - Z 2 )1>) 



In this case the action (I4.39[) becomes 5 complex matrix model instead 
eight real matrices Ya^, Yb^ and one complex matrix ip. This form may seem 
cumbersome but it still more comfortable for numerical simulations. The next 
step is to define the estimator for the average values of interest and to develop 
some numerical parameters in order to analyze the numerical results. 



4.2 Definition of the observables 

Calling (if), Zi) a configuration of the fields (where 2 = 1,2 for the two dimen- 
sional case and i — 1, • • • ,4 for the four dimensional one), the probability to 
encounter this configuration is given by 

e -s[y>,Zi)) 

P[( 1 j,,Z i )] = — (4.41) 

where S[(i/), Zi)] is the action (14.391) or (14.401) of the system evaluated in the 
particular configuration and Z is the partition function: 

Z = D[(^,Zi)]e- s[ ^' z ^ (4.42) 

D[(ip, Z^] denotes the integration over all field configurations. The expectation 
value of the observable O is defined by the expression: 

(O) = / DM, Zi)] e - (4.43) 



70 



Chapter 4. Preparing for numerics 



Following Monte Carlo methods, will be produced a sequence of configurations 
{(ip, Zi)j}, j = 1,2, ••• , Tmc an d evaluated the average of the observables 
over that set of configurations. The sequences of configurations obtained, a 
Monte Carlo chain, are representations of the configuration space at the given 
parameters. In this frame the expectation value is approximated as 

1 Tmc 

Oj (4.44) 
i=1 

where Oj is the value of the observable O evaluated in the j-sampled configu- 
ration, (ip, Zi)j, Oi = 0[(if), Zijj}. 




4.2.1 Energy and specific heat 

The internal energy is defined as: 

E(Q,fi,a) = (S) (4.45) 
and the specific heat takes the form 

C(Q, f i,a) = (S 2 )-(S) 2 (4.46) 
These terms correspond to the usual definitions for energy 

£(fi,/i,a) = --— (4.47) 

and specific heat 

C(n,^,a) = — (4.48) 

where Z is the partition function. It is very useful to compute separately the 
average values of the four contributions: 

Sjr^Zi) = Tr Cp (4.49) 

SvMZi) = Tr A/ (4.50) 

S Vl (ij),Zi) = Tr£ yi (4.51) 

S D {^,Zi) = Tr(C Dj £ Dj ) (4.52) 

Where i,j — 1, 2 or i, j — 1, • • • , 4. The corresponding expectation values to 
(gagD - fl4~52j) are 

F(Q,fi,a) = (S F ) (4.53) 

V (n,ti,a) = (S Vl ) (4.54) 

V^^a) = (S Vl ) (4.55) 

D(n,n,a) = (S D ) (4.56) 

V = Vo + Vi (4.57) 
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4.2.2 Order parameters 

The previous quantities are not enough if we want to measure the various 
contributions of different modes of the fields to the configuration tp, Z^. There- 
fore we need some control parameters usually called order parameters. As a 
first idea we can think about a quantity related to the norms of the fields 
for example the sums ^2 nm |Vw| 2 , ^2 nm \ Zinm\ 2 of all squared entries of our 
matrices, this quantity is called the full power of the field [HD E2] and it can 
be computed as the trace of the square: 

Vl = Tr(H 2 ) (4.58) 
Zl = Tr(|^| 2 ) (4.59) 

This parameter turns to be very useful to find the region where the con- 
figurations are considered disordered for which we expect to have a random 
distribution of the modes ((pi) ~ 0, (Zf a ) ~ and the so called ordered phase 
in which we have ((pi) » 0, (Zf a ) » 0. The full power of the field can 
be used to determine for which independent parameters there is a transition 
between the disordered and ordered regime. However (ip a ) alone is not a real 
order parameter because does not distinguish contributions from the different 
modes but we can use it as a reference to define the quantities: 

N 

,r? l„ |2 



£ 

n=0 
N 

Zl = $> mn | 2 (4.60) 

ra=0 

Referring to the base ( I4.22p it is easy to see that such parameters are connected 
to the pure spherical contribution. This quantity will be used to analyze 
the spherical contribution to the full power of the field. We can define the 
corresponding susceptibility as: 

X V0 = (<Po) ~ (Vo) 2 

X z M = (Zl)-(Z l0 ) 2 (4.61) 

We can generalize the previous quantity defining some parameters ipi in such 
a way they form a decomposition of the full power of the fields. 

= + E Z l = Zl + Yl Z l (4.62) 
i i 

Following this prescription the other quantity for I > can be defined as: 

i i 

= 7", Km(l - 8 nm )\ 2 , Z\ = 22 \ z lnm0- ~ 8nm)\ 2 (4.63) 
n,m=0 n,m=0 
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If the contribution is dominated from the spherical symmetric parameter we 
expect to have ~ (v^o); i^ia) ~ (^fo)- m ^his regime for the uniform 

ordered region we have (</?q) ~ 0, (Zf ) ~ and in the disordered phase we 
expect ((pf) ~ 0, (Zf t ) ~ for all I. With (I4.63[) we can define the order 
parameters ipf, Zf ± as a particular case I = 1 we have 

<fl = Ko| 2 + Ki| 2 , Z\ = \z il0 \ 2 + \z m \ 2 (4.64) 
and its susceptibility, 

X*i = - (^i) 2 , Xz a = {ZD - (Za) 2 (4.65) 

In the next simulations we be evaluated the quantities related to I = and to 
/ = 1 as representative of those contribution where the rotational symmetry 
is broken. A region characterized by (<pf) » 0, (Z^) » is called [SI] 
non-uniform, such region feature a large contributions in which the spherical 
symmetry is broken. Using higher / in (14. 63p we can analyze the contributions 
of the remaining modes, but turns out that the measurements of the first two 
modes are enough to characterize the behavior of the system. 
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Numerical results 



This last chapter is devoted to present the numerical results of the 4- dimensional 
and 2-dimensional model constructed and approximated in the previous chap- 
ters. In the first part will be discussed the behaviors of the defined observables 
of the full 4- dimensional model varying the parameters. Meanwhile in the sec- 
ond part will be showed the result for the 2-dimensional model for the same 
values of parameters and compared with the full model results. 



5.1 Four dimensional case 

Now we are ready to discuss the results of the Monte Carlo simulation on 
the approximated spectral model described in the chapter 2 and discretized in 
the previous chapter. As a first approach to Monte Carlo simulation of such 
kind of model, we use some restrictions on the parameters. Will be used the 
spectral action around its minimum (I4.40p in order to simplify the calculation, 
in this frame in the action is symmetric under the transformation — > — // 
so will be used the condition /z > and /i 2 > 0. Will be explored the range 
/i G [0, 3], this interval was chosen in such a way to show a particular behavior 
of the system for fixed Q. The parameter Q appears only with its square and is 
defined as a real parameter, therefore for the Q too we require Q > 0. Beside, 
if we refer to the scalar model, is possible to prove [H] using the Langmann- 
Szabo duality [T7], that the model can be fully described varying Q in the range 
[0, 1], for higher fl the system can be remap inside the previous interval. In the 
present model the L-S duality does not hold any more, but forced by limited 
resource, we conjecture that the interval Q e [0, 1] still enough to describe the 
system. The last parameter to consider is a, it is connected to the choice of 
the vacuum state, of course, the range of a is [0, 2ir]. The study of the system 
varying this parameter is quite important; from a theoretical point of view, 
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because is related to the vacuum invariance and for numerical reasons too. In 
fact, in the action of both 4-dimensional and 2-dimensional model appear some 
contributions ~ (sina)/f2 which seem to diverge for Q = 0, numerically we 
have verified that is an eliminable divergence and the curves of the observables 
can be extended in zero by continuity. Nevertheless, this factor slow down the 
computations due to its high value close to Q — 0. Studying the dependence on 
a we can conclude that in the limit N — > oo the observables are independent 
from a, therefore for our purposes a will be fixed equal to zero avoiding the 
annoying terms. In general for each observable are computed the graphs for 
iV =5, 10, 15, 20 matrix size. 

5.1.1 Varying a 

We start looking at the variation of energy density and the full power of 
the fields density for fixed fi and Q, varying a £ [0, 27r]. As representative 
here will be presented the graphs for \i — 1, Q — {1,0.5} but we obtain the 
same behaviors for any other choice of the parameters allowed in the range 
considered. All tree graphs show an oscillating behavior of the values, this 
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Figure 5.1: Total energy density and full power of the fields density for <p^, Z"^ a (from the 
left to the right) fixing fi = 1, fi = 1, varying a and N. N — 5 (blue), N — 10 (purple), 
N = 15 (brown), N = 20 (green). 

oscillation is present in all other quantities measured. The amplitude of this 
oscillation becomes smaller and smaller increasing the size of the matrix and 
this is true for all the quantity measured. The same trend is described in 
fig l5.2l in which position of the maximum are different but the amplitudes 
become smaller increasing TV 

This results allow us to consider a = for all next graphs, since we are 
interested in the behavior of the system for iV — > oo. This occurrence simplify 
all the next simulations thanks to the vanishing the of terms ~ (sina)/f2 
appearing in the discretized action. Beside, such results induce us to reckon 
the parameter a as connected to the remaining invariance of the vacuum state 
for the exact model. 
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Figure 5.2: Total energy density and full power of the fields density for ip"^, Z^ a (from the 
left to the right) fixing /i = 1, Q. = 0.5, varying a and N. 

5.1.2 Varying Q 

Now we will analyze three cases in which /i is fixed to 0,1,3, a is zero and 
we vary Q e [0, 1]. Will be showed the results of simulation for matrix size 
of 5, 10, 15, 20. In the rest of this chapter we ignore for the computation of 
E,D,V,F the prefactor (1 + fi 2 ). In this way we focus our attention to the 
integral as the source of possible phase transitions. The graphs in fig J5.3l show 




Figure 5.3: Total energy density and the various contributions for fi = 1. From the left 
to the right the density for: E, V, D, F. 

the total energy density (S)/N 2 and the various contributions: the potential 
V/N 2 , the Yang-Mills part F/N 2 and the covariant derivative part D/N 2 , for 
jj, — 1 a. — defined in the previous chapter. There is no evident discontinuity 
or peak, apart a small deviation visible in the origin, but we can interpret it 
as a finite volume effect. Increasing the size of the matrices the curves and the 
apparent discontinuity close to zero become smoother fig J5.4l 
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Figure 5.4: Total energy density and the various contributions for n = 1, a = varying f2 
and N. From the left to the right E, V, D, F with N = 5 (blue), N = 10 (purple), iV = 15 
(brown), N = 20 (green). 

Comparing the energy density and the various contributions fig J5.5l we no- 
tice that the contributions between F and V balance each other and the total 
energy follows the slope of D, this behavior continues increasing the size of 
the matrices fig J5.5l 




Figure 5.5: Comparison of the total energy density and the various contributions for 
H = 1, a = 0. E (blue), F (brown), D (green), V (purple). With N = 5 (left) and N = 20 
(right). 
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Figure 5.6: Specific heat for (i = 1. 

The specific heat density defined as ((S) 2 — (S 2 ))/N 2 shows fig J5.6l a peak 
in Q — 1, this peak increases as N increase. This behavior is typical of a phase 
transition and the peak is not clear for small N due to the finite volume effect. 
In order to gain some informations on the composition of the fields we look at 
the order parameters defined in the previous chapter. Starting from ip field, 
in the figure EE it is showed the graphs for ip 2 , ipl and <p\ for N = 5. 
The three values ip 2 , <pl and (p 2 seem essentially constant, comparing the three 




Figure 5.7: From the left to the right ip 2 a , y\ and ip\ for ^ = 1 varying ft 

graphs fig J5.7l it is easy to see the dominance of the spherical contribution ^ 
to the full power of the field. 
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Figure 5.8: On the left comparison of ip\ (blue), tp\ (purple) and <p\ (green) density. On 
the right comparison of Z^ a (blue), Zq (purple) and Zq X (green) density. 



The behavior of the Zq fields is different, referring to figure [5781 where are 
showed the quantity Z^ a , Zq and Zq X over N 2 , the spherical contribution 
becomes dominant approaching to Q = starting from a zone in which the 
contribution of Z^ 
only the graphs for Zq g , Zf 71 
the other related graphs appear compatible to the Z case. The dependence 
of the previous quantities on N are showed in the following graphs fig J5.9l The 



00 and Z\ x are comparable. For brevity will not be showed 
0a , ^ 00 and Zq X but taking in account the statistical errors 





Figure 5.9: Starting from the up left corner and from the left to the right the densities 
for <p*, <fo, (fil, ZQ a , Zq Q and Z^ x for fj, = 1 varying f2 and N. 

values of the quantity of all the previous parameters decreases with iV but the 
dominance of the (po on the total power of the field is independent by N. The 
peak related to Z decrease, but if look at the single graph for N = 20 the 
spherical contribution approaching the point Q = features a peak. 

Now we will analyze the model for /i = 0; fig J5. 101 shows the graphs for total 
energy density and the contributions V, D, F. The slope of the total energy 
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density seems to be constant. The D contribution and the F do not balance 
each other like in the previous case, but all the three contributions balance 
among them self to produce a constant sum. The specific heat density shows 




Figure 5.10: Total energy density, various contributions and the comparison among them 
for fi — varying il and N. From the left to the right E, V, D, F and comparison. 
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Figure 5.11: Specific heat density for /j, = varying £1 and N. 

fig J5.11l again the peak in Q = 1 as N increase. For the other quantities tp^, (fl, 
ip\ and Zl a , Zq , Zq X we have the same behavior fig J5.14l of \i = 1 case, except 
for the oscillation appearing in the Z^ a , Z^ graphs close to zero, anyway it 
appears only for N = 5. 
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Figure 5.12: Starting from the up left corner and from the left to the right the densities 
for <p*, <fo, <Pi, ZQ a , Zq Q and Z^ x for fi = varying f2 and N. 

A complete different response of the system is described in the graphs for 
\i = 3, as we can see from fig J5.13l The slope of total energy density is very 
similar to the F component instead D. Beside, appears a maximum around 
Q = 0.4 for iV — > oo. This dramatic change of the graphs might be interpreted 
as consequence of a phase transition ruled by the parameter /i, actually in the 
next section we will find a peak in the specific heat density for some fixed Q 
and varying \i e [0, 3]. 



□ l.L' 




Figure 5.13: Total energy density, various contributions and the comparison among them 
for /i = 3 varying and N. From the left to the right E, V, D, F and comparison. 
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Specific heat density displays fig J5.14l this change too, in fact instead the 
peak in O = 1, it appears in the opposite side of the studied interval in fl = 0. 
This peak too, due to its grows increasing N indicates a phase transition. 




Figure 5.14: Specific heat density for /z = 3 varying f2 and N. 




Figure 5.15: Starting from the up left corner and from the left to the right the densities 
for <p*, <fo, (fll, ZQ a , Zq and Z^ x for fi = 3 varying f2 and N. 

The fig J5. 151 describes the behavior of the order parameters densities (p^, ip^, Lp\ 
and Zq q , Zq , Zq 1; they have a similar aspect to the previous relative graphs. 
For the ip field the spherical contribution remains dominant, beside in the 
graph appears a deviation from the constant slope this deviation is evident 
for N = 5 but still there for higher N. The order parameters for Z display a 
peak close to the origin without oscillations even for N = 5. This maximum 
for higher iV does not move closer to the origin, in other words, this shift in 
not due to the finite volume effect. Even for Zq X graph appears a deviation 
from the constant slope, a small peak which becomes shifted and smoother for 
higher N 
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5.1.3 Varying fi 

In this section is analyzed the response of the system varying fi 6 [0, 3] where Q 
is fixed to 0, 0.5, 1 and a is always zero. We start displaying the graphs fig J5.16l 
of the total energy density and of various contributions for Q = 0. There is no 
evident discontinuity but appears a peak in the total energy density around 
jj « 2.5 for iV = 20. Comparing all the contributions is easy to notice that 
the slope of the total energy is dictated by the curve V of the potential part. 




Figure 5.16: The total energy density and the various contributions for f2 = varying fj, 
and TV. From the left to the right E, V, D, F an comparison with N = 5 (blue), N = 10 
(purple), N — 15 (brown), N = 20 (green). For the comparison: E (blue), V (purple), D 
(brown), F (green). 




Figure 5.17: Specific heat density for ft = varying /i and N. 

As mentioned before, the specific heat density fig J5.17l features a peak 
around \i ~ 2.5 for iV = 20 and again, due to this behavior as N increase, we 
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could relate this peak to a phase transition. The plots for the quantities y>\ 
and ipl denote a strong dependence on fi, in particular the slope of tp\ seems 
mostly linear, for ip\ we notice a similar behavior but the slope is no longer 
linear. Comparing the three graphs fig l5.18l we deduce that close to the origin 
the non spherical contribution is bigger the spherical one, increasing \i this 
situation capsizes and y?g becomes dominant respect ^>\. The behavior of the 




Figure 5.18: Starting from the up left corner and from the left to the right the densities 
for <£q, ipl, Zq q , Zq Q and Z^ x for £1 = varying \i and N. 

Zq fields is quite different, referring to figure 15.181 the spherical contribution 
is always dominant for the all interval \i G [0,3]. The curves for Z^ a , Zq Q are 
compatible to the constant slope, for Zq X we have the same dependence on 
fi in particular there is a smooth descending step, however this step becomes 
smoother for bigger N. It is behooves to say that due to some cancellations 
effects the statistical errors are quite big and they can hide some dependence, 
anyway this results tell us about the dependence of the order parameter for 
Zi and in general of the system, on the two choice fi = or fi ^ 0. 

Now we will analyze the model for fi = 0.5; as fig l5.19l shows the graphs 
have a different slope comparing to the previous case, the maximum of total 
energy density follow the one of the V component. If we focus ourself only on 
the total energy graph and we compare it with the one for = 0, we notice 
a shift of the maximum for each N. In particular in fig J5.19l some maximum 
are moved outside the studied interval. We can find this shift very clearly 
looking at specific heat density graph fig J5.20} we find again the peak as N 
increase but it is shifted around /i ~ 3.3. The graphs fig l5.23l for have 
the same behavior of Q = 0.5 case, excluding some fluctuations close to the 
origin for N = 5 due to the finite volume effect, ipf graph displays an almost 
constant curve. However, close to the origin, the spherical contribution and 
the first non spherical one are comparable. The introduction of fi 7^ creates, 
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Figure 5.20: Specific heat density for = 0.5 varying /i and TV. 



in the Z fields order parameters fig J5.214 a dependence similar to the graphs 
for the if); the full power of the field density and the spherical contribution are 
no more constant and they grow increasing /x. Even in this case the spherical 
contribution is always dominant excluding the region around fi = 0. 
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The last set of graphs for the 4- dimensional model are obtained fixing 
Q — 1, due to the vanishing of pref actor in front of the Yang-Mills part of 
the action the F contribution is always zero. The following diagrams for the 
energy and contributions show the absence of the previous peak and comparing 
again them with former graphs they seem a sort dilatation. The specific heat 




Figure 5.22: Total energy density and the various contributions for f2 = 1 varying /i and 
N. From the left to the right E, V, D. 

density does not show the peak in zero any more fig J5.23l and the curves does 
not show any particular point as N increase, actually the peak can be found 
for higher /x. At last in fig J5.24l we found a behavior of the order parameters 
density for Zq and i/j fields similar to the former graphs for Q = 0.5 and they 
are compatible with a dilatation of the previous diagrams. 
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Figure 5.23: Specific heat density for ft = 1 varying fi and N. 




Figure 5.24: Starting from the up left corner and from the left to the right the densities 
for <fl, <Po, <Pi, z la-, z oo and z oi for SI = 1 varying fi and N. 

5.2 Two dimensional case 

The main aim of this chapter was to show the numerical result for the 4- 
dimensional model, beside it is very interesting to analyze the 2-dimensional 
model (which was used as a first test of the program) described by the dis- 
cretized action (I4.39P and note the differences with the full model. Will be 
used the same set of parameters used for the full model and due to the same 
dependence on a all the simulation will be conducted with a = 0. 

5.2.1 Varying Q 

The behavior of the energy density for \x = {1, 0} is the same of the full model, 
in the graphs fig 15.251 are plotted the total energy density and the comparison 
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of the various contributions for /x = {1, 0}. As we can see there is no evident 
difference from the 4- dimensional case; even in this case the total energy for 
/x = 1 follows the slope of the D contribution and the total energy for /x = 
is constant. Even the various contributions have the same trend of the full 
model; for /x = 1 the V and F have an "opposite" slope and for /x = they 
balance each other in order to obtain a constant sum. The specific heat density 
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Figure 5.25: Total energy density and contributions for fi = 1 (left side) and \i = (right 
side) varying Q and N. Where E (blue), V (purple), D (brown), F (green). 



shows fig |5.26l some differences: for fi = the peak is no more located around 
Q = 1 but it appears, not so clearly, around zero and the specific heat for 
H = does not show any particular point. In the figure fig J5.27l are showed 




Figure 5.26: Specific heat densities for fj, = 1 (left) and \x = (right) varying tt and N. 

the graphs for ip 2 a1 y>\ and tp\ and Z^ a , Zq Q and Zq X for /x = 1. For brevity 
are showed the graphs only for Zq, however even for the 2-dim model taking 
in account the statistical errors they are compatible with the Z\ field related 
graphs. In this results there are no significant changes with the 4-dim model, 
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Figure 5.27: Starting from the up left corner and from the left to the right the densities 
for (pi, <fo, ZQ a , Zq Q and Z^ x for = 0.5 varying f2 and N. 

indeed comparing the three graphs fig J5.27l it easy to see the dominance of the 
spherical contribution to the full power of the field and for the Z Q fields the 
spherical contribution becomes dominant approaching to Q = starting from 
a zone in which the contributions of Zq and are comparable. Beside, we 
notice a decrease of values of the quantiles related to if) with N, nevertheless 
the dominance of the ipo on the total power of the field remains independently 
by N. The peak appearing in Z^ a and Zq decrease, but if look at the single 
graph for fixed N the spherical contribution approaching the point £1 = still 
dominant. The graphs of the previous quantities for /x = show the same 
slope as the full model so will be not presented here. 

A radical change emerges for the case \x = 3 as in the full model, but the 
plots of the energy density are very different fig |5.28l From the comparison of 
the various contributions with the total energy density for iV = 20, we deduce 
that the total energy density seems to follow the D contribution instead F as 
in the 4-dim model. This can be justified by the presence of additional terms 
in the full model coming from the Y-M part. 
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Fi gure 5.28: Total energy density and contributions for \x = 1 varying Q and N. From 
the left to the right E, V, D, F and comparison. 
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Figure 5.29: Specific heat density for /z = 3 varying fi and N. 

In this case the specific heat density does not shows fig l5.29l the peak for 
Q — 1, but we discover a step around \i ~ 0.5 this seems to becomes sharper as 
N increase. Regarding the order parameters ip 2 a , (p^, <pf and Z^ a Zq , they 
have similar behavior of the full model with some more evident oscillations 
around the origin. 



5.2.2 Varying fi 

The resemblance between the two models is more evident when we compare 
the plots obtained fixing Q = {0, 0.5, 1} and varying \i. The graphs in fig J5.31l 
show the total energy density and the various contributions for Q = {0, 0.5, 1}. 
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Figure 5.30: Starting from the up left corner and from the left to the right the densities 
for Vl, Vo, Vi, z la-, z oo and z oi for A* = 3 varying fl and N. 

It is very clear the similarity with the full model; for Q = appears a peak in 
the total energy density around /i ~ 2.3 for iV = 20, for Q = 0.5 and Q = 1 
the peaks are shifted. The various contributions follows the same behavior of 
the full model where, in the case $1 = 1, the F component is absent due to the 
vanishing of the prefactor in front of the Y-M part of the action. 
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Figure 5.31: Total energy density and contributions for fi = I, a = varying fl and N. 
From the left to the right E, V, D, F and comparison with N = 5 (blue), N — 10 (purple), 
N = 15 (brown), N = 20 (green). 
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Figure 5.32: Specific heat density for O = (left) and O = 0.5 (center), = 1 (right) 
varying \i and N. 



The resemblance persists in the specific heat density too, in fig J5.32l the 
peaks are located as the 4-dim model and for Q = 1 does not appear any peak. 
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Fi gure 5.33: From the left to the right the densities for Z^ al Zqq and Zq^ for = varying 
H and N. 

At last in fig J5.33l is displayed the behavior of the order parameters density 
of Z for Q = 0, even in the 2-dimensional case the slopes of Z^ a , Zq densities 
are constant within the errors bars. The other orders parameters are identical 
to the full model related graphs and are compatible with a dilatation of the 
previous diagrams for Q = 0.5. 
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In this thesis we have presented a study of a spectral action model con- 
structed in order to extend to the Yang-Mills theories the non- commutative 
Wulkenhaar-Grosse model. The main aim of this work was to test a first Monte 
Carlo approach based on a non-perturbative regularization method alternative 
to the standard lattice regularization, which in our case is not suitable due to 
the oscillator factor of the Moyal product. We have performed Monte Carlo 
simulations and obtained the values of the defined observables varying the pa- 
rameters of the system. Despite the complexity of the approximated spectral 
action considered here we were able to obtain some reliable numerical results, 
we can conclude that a numerical approach of this kind of model using the ma- 
trix Moyal base seems feasible. The specific heat density shows various peaks 
indicating phase transitions, in particular studying the behaviors for some fixed 
/i we found a peak around Q = for ji = {0, 1} and peak in = 1 for /i — 3, 
beside we notice a huge change in the energy density and in its contributions 
between the cases [i = {0, 1} and /i = 3. Others peaks in the specific heat 
density can be found varying [i and fixing Q, the graphs show that increasing 
Q the peak in the specific heat, starting from /i ps 2.4 for Q = 0, is moved 
towards higher \i. The order parameters introduced show a strong dependence 
on the occurrence of f2 = or f2 ^ 0. Referring to the fixed /i graphs we found 
a peak in the spherical contribution for the gauge fields Zj, we can interpreted 
this slope as a sort of symmetry breaking introduced by Q ^ 0. Additionally, 
varying fi and fixing Q the other parameters display an increasing slope with 
fi for all fields and all situation but one; the graphs of the order parameters 
concerning Z 0a , Z Q0 for Q = show a constant behavior. As an aside result 
we have ran the same computations on the simpler 2-dimensional model, from 
the comparison despite many similarities in the results varying fi, we note a 
complete different trend of the specific heat varying Q, the peaks present in 
the full model does not appear any more in the 2-dimensions one. From this 
comparison we deduce that the dimension is a crucial parameter of the model. 
About the prospective, there is a lot of work to do; the natural next steps in 
the numerical study of this model, could be the computation of the transition 
curves in order to separate the phase regions and classify them using eventually 
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some additional order parameters. Our treatment, forced by limited resource, 
was conducted conjecturing that the system can fully described varying Q in 
the range [0, 1] but since the L-S duality does not hold any more in our case, 
will be very useful to extend this range. Actually, the computed graphs does 
not show any periodicity in Q G [0, 1] so there are some hints to infer that 
the range [0, 1] is not enough, however only a direct computation will clarify 
this point. At last, in analogy with the matrix model will be very interest- 
ing to do not require any more the condition fj 2 > 0, in order to implement 
this change we have to conduct the calculation no more around the minimum 
of the actions. The problems that can occur are related the thermalization 
process, which can be much more problematic due to the eventual emerging 
of local minimum. The expansion of the parameters space, together with the 
classification of the different phase regions allows us to compare our model 
with the results of the simulation conducted on the fuzzy spaces, looking in 
particular about the occurrence of the so called non-uniformly ordered phase 
which is connected with the UV/IR mixing. Since we have constructed our 
model starting from a renormalizable one this study is very desirable. 
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Appendix A 

Axioms non-compact geometries 



In this appendix will be introduced the definition of non- compact spectral triple 
and the modified set of axioms required for such spectral triple. In the last part 
of the chapter will be hinted a more rigorous treatment of the spectral triple 
used in chapter 3, extending it in the setting of non-compact triple. 

A.l General remarks 

A compact spectral triple is a collection (A, V, %) where A is a unital pre-C*- 
algebra, % is a Hilbert space on which A is represented by bounded operators 
and D is a Dirac selfadjoint operator on A, with resolvent -Rd(A) = (V — A) -1 
compact and the commutator [D, a] is bounded on % for every a £ A. In 
the non-compact case this definition is modified to allow A to be no longer 
unital. The usual way to construct a non-commutative spin manifolds, for the 
unital case, is to start from an operatorial version of ordinary spin geometry 
that can be generalized to a non- commutative manifold. This generalization is 
obtained using a reconstruction theorem in order to recover all of the structure 
from the abstract geometry over a suitable algebra [28| l35| 168] . However, there 
is no proper direct reconstruction theorem for non-compact Riemanian spin 
manifolds, so the question of axioms for non-unital spectral triples is guided 
by the resemblance to the compact case, requiring a collection of axioms for 
not necessarily compact non-commutative manifolds. Of course, such list of 
conditions should be compatible with the previous axiomatic framework and be 
fulfilled by non-compact commutative manifolds. Modified systems of axioms 
for non-compact spectral triples has been formulated by A. Rennie [69] for 
the commutative case and by Gayral et al. [55] in a fully non-commutative 
framework on the Moyal plane. A very important aspect for these formulations 
is the introduction of various different types of non-commutative dimension 
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[70] . Indeed, the non-compact spectral triple analyzed will be characterized 
by different metric and KO-dimensions. When A is not unital, we need to 
choose a preferred unitization B. For compact or unital spectral triples it 
is further required that the operator T> has a compact resolvent, or in other 
words, (T> — A) -1 must belong to K(~H) for A ^ spV, where K(7i) denote the 
compact operators on "H. Consequently T> must have discrete spectrum of finite 
multiplicity. Since this is clearly not the case of Dirac operator Tp on M. d , in 
the nonunital case we only demand (35] that a(T> — A)^ 1 is compact for a £ A. 
This condition ensures that the spectral triple (A, V, H) corresponds to a well- 
defined (^-homology class [73] for an, in general, non-compact non-commutative 
geometry. Now we analyze some modified conditions which spectral triples are 
required to satisfy. We use a regularity condition (FTJ in terms of the operator 
(V) = (V 2 + 1)2, in this way we avoid the hypothetical zero modes of the 
Dirac operator. Can be shown that the previous regularity condition coincides 
with the original definition of Connes and Moscovici [3S]. The condition of 
finiteness, in the unital case, is that the smooth domain 1-L°° of T> in % can 
be finitely generated by a projective left module over the unital algebra A, or 
"H°° ~ A m p for some projector p in M m (A) with a suitable m. In the case 
of M. d , using the commutative or the Moyal product, the module of smooth 
spinors is free, since the triviality of the spinor bundle. 

When A is nonunital, in order to obtain a projective ^4-module, we have 
to find the projector p in a matrix algebra over the preferred compactification 
B, therefore this condition for the nonunital case should then require that T-L°° 
densely contains a pullback of a finite projective i3-module to A. The reality 
axiom is the same as the compact case, thus is the existence of an antilinear 
conjugation operator J on "H such that a — > Ja* J" 1 , which define an another 
representation of A on "H commuting with the previous one. The required 
algebraic properties of the real structure J are the same as the compact case. 
Even for the first order condition, there is no need to modify the axiom in the 
non-compact case, we have: 

[[£>, a], JtfJ- 1 } = 0, for all a, 6 e A. (A.l) 

It is convenient formulate the orientability condition in the compact case as 
the property that the spectral triple (A, H,T>) contains an algebraic version 
of a volume <i-form. This algebraic version consists of an Hochschild <i-cycle 
c, which is a sum of terms of the form (do <S> b^) <S> «i <E> ■ ■ ■ <8> ad satisfying the 
cycle property be = 0. Defining 6q as an element of the opposite algebra A" 
with bo € A, we can represent the Hochschild <i-cycle by bounded operators: 

^((ao ® fig) ® ai ® ■ • • ® a d ) = a J6J J _1 [X>, a t ] ■ ■ ■ [V, a d }. (A.2) 

on which we impose the orientation 7T£>(c) = X-, where x is the given Z 2 -grading 
operator on %. We just use x = 1 if <i is odd and in the even case, for ordinary 
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spmors, we use 

X = (-0 m 7 1 7 2 ---7 am (A.3) 

For nonunital algebras, we ask the same condition taking the cycle over the 
unitization B rather than A itself. 



A. 2 Modified conditions for nonunital spec- 
tral triples 

By a real non-compact spectral triple of dimension k we connote the data 
(A, B , "H, D, J, x), where A is an algebra, in general nonunital, acting faithfully 
via a representation denoted by n on the Hilbert space H, B is a preferred 
unitization of A acting on the same Hilbert space and D is an unbounded 
selfadjoint operator on % such that [D, a] for each a in B extends to a bounded 
operator on H. In addition, J and x are respectively an antiunitary and a 
selfadjoint operator, such that x = 1 when k is odd and otherwise x 2 = 1> 
Xa = ax for a G A and T>x = ~xP- This objects are asked to satisfy the 
conditions below. 



0. Regularity 

The bounded operator , n(a) and [D, n(a)] for all a G B belongs to H^Li Dom(5 n ), 
with 5(T) = [1(D), T] and (D) = (D 2 + l)s 



1. Compactness 

The operator 7r(a)("D — A) -1 is compact for all a E A and A G sp(D), where 
sp(D) is the spectrum of D. 



2. Spectral dimension 



For any element b of the algebra ^q(A) generated by 5 n (7r(A)) and 5 n ([D, tt(A)]), 
the function (b{ z ) — Tr(6("D) _z ) is well defined, holomorphic for %K{z) large 
and analytically continues to C/Sd, where Sd C C is a discrete set called di- 
mension spectrum. The dimension spectrum is said simple if all the poles are 
simple, finite if there is a k G N such that all the poles are of the order at most 
k and if it is finite. 
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3. Metric dimension 

The metric dimension is defined as d — sup{9i(z), z G Sd}, the operator 
7r(a){V)~ d belongs to the Dixmier class C^H) for any a £ A. For any Dixmier 
trace the map a G A+ — > Tr u ,(Tr(a)('D)~ d ) in non vanishing. 

4. Finiteness 

The algebra A and its preferred unitization B are pre-C*-algebras. Each one 
is a *-subalgebra of some C*-algebra, stable under holomorphic functional 
calculus. The space of smooth spinors 

oo 

U°° = p) n n , where U n = Dom{V n ) (A.4) 

n=0 

completed with the norm | |£] | n = | |£| | 2 + 1 \T> n ^\ | 2 , is a finitely generated projec- 
tive A- module pA m for some m G N and some projector p = p 2 = p* G M m (B). 
The scalar product (•, •) on "H°° can be obtained composing the Dixmier trace 
and the induced structure (•, ■) : W° x W° -> A 

Tr UJ ((^v}(V)- d ) = (Z,v), C^e^ 00 (A.5) 

5. Reality 

There is an antiunitary operator J on T-L, such that: 

[n(a), Jn{b*)J- 1 } = 0, for all a,b G B (A.6) 

b — )> Jb* J~ l is a commuting representation on H of the opposite algebra A". 
Beside, the operator is required to satisfy J = e and JV = e'VJ and also 
for the even case Jx = t"xJ , where e, e', e" G { — 1, 1} given as a function of 
k Q G Z 8 . This function is represented by the table: 



ko mod 8 





1 


2 


3 


4 


5 


6 


7 


e 


1 


1 


-1 


-1 


-1 


-1 


1 


1 


e' 


1 


-1 


1 


1 


1 


-1 


1 


1 


e" 


1 




-1 




1 




-1 





(A.7) 



This ko is usually referred as the KO dimension of the spectral triple. 
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6. First order 

The bounded operators [X>, 7r(a)] also commute with the opposite algebra rep- 
resentation: 



7. Orientation 

There is an Hochschild d-cycle c on B, with values in B ® B'°. Such <i-cycle 
is a finite sum of terms like (a <g> b°) <S> ai <8> • • • <8> Qfc, which can be naturally 
represented by operators on % using 7Tx>(c) in formula (lA.ip . The volume form 
7Td(c) must solve the equations: 



A geometry is called connected or irreducible if the only operators commuting 
with A and T> are the scalars. 

Poincare duality, for a non-compact orientable manifold M, is usually ex- 
pressed as the isomorphism between the compactly supported de Rham coho- 
mology and the homology of M mediated by the fundamental class [M]. We 
shall actually leave aside the final condition of Poincare duality, since it is not 
central in the present work. 

A. 3 Non-compact harmonic Moyal case 

Referring to the spectral triple (A4, V4, H4) defined in the chapter 3 it is pos- 
sible to extend the previous triple to a non-compact spectral triple of metric 
dimension 4 and KO dimension 0. In this section will be just listed the theo- 
rems about compactness, boundness, regularity and dimension remanding to 
[57] for the complete proofs. 

Compactness and Boundness 

As first step we define a preferred unitisation £>4 of A4 as the space of smooth 
bounded function on IR 4 with all partial derivatives bounded. The Moyal 
product extend to #4 and A± C B4 is not a dense two-side essential ideal [55] . 
The algebra £> 4 contains the plane waves e %kx and constant functions but no 



[[£>, tt(o)], Jir(b*)J ] = 0, for all a, b e B. 



(A.8) 



7Td(c) = X even case 
7T£>(c) = 1 odd case 



(A.9) 
(A.10) 
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other non-constant polynomials. Can be proved [55] that £> 4 is contained in 
the C*-algebra: 

A 4 = {TE S'(R 4 ) : T*f G L 2 (R 4 ) for all / G L 2 (R 4 )} (A.ll) 

completed using operator norm. Defining the action n = L* of £> 4 on "H 4 as 
the componentwise left star multiplication: 

*:A 4 xH^H 4 , (/,V)^^*(/)V = /*^ (A. 12) 

For / G B4 the operator 7r(/) is bounded. The commutator [1*4, 7r(/)] is 
computed using the relation (I2.25P 

V A {fk^)-fk{V^) = (iV + nr^drf))^ (A.13) 

with / 6 B4 and t/> G Dom(£> 4 ) D Dom(£> 4 )7r(/). We deduce that for <9 M / G 
fi 4 the commutator is extended to a bounded operator. The compactness is 
verified taking in account that [V>\ + 1) _1 is already a compact operator on 
H4, as bounded operator K is compact if and only if KK* is compact too. 
In particular, we have that the operator (T> 4 + i)' 1 is compact, beside the 
definition resolvent equation implies that (X> 4 + A) -1 is compact for any A in 
the resolvent set of X> 4 . Whence 7r(/)(£> 4 + A) -1 is compact for / G £> 4 . 

Regularity and dimension 

As in the compact case it is useful for the non-compact too, to define the 
unbounded operators R and L as 

R(T) = [VlT](V 4 )-\ L(T) = {V A )- l [Vl,T\ (A.14) 

Can be proved that they satisfy f]^Li Dom5 n = f]^ n=1 DomR n L m . Using this 
derivation can be demonstrated [57] the following proposition: 

Proposition 1. Both 7r(/) and [D4, 7r(/)] belong for any f G B4 to f]^ =1 Dom5 n , 
where 5 n are defined by the axiom 0. 

For the dimension spectrum we have the following theorem: 

Theorem 1. The spectral triple {A±,T '4, ^4) has dimension spectrum Sd= 
4 — N and therefore metric dimension 4- 
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Error estimation and update 
algorithm 



In this appendix will be introduced the treatment of the statistical error in 
Monte Carlo simulation and will be discussed three methods to estimate the 
statistical errors. In the end of the section will be discussed some issues con- 
cerning the estimation of the autocorrelation time and the implemented update 
algorithm. 

We start defining the mean O over a subset of {OijJ^f of T elements as 
a simple arithmetic mean over the Markov chain: 



where T MC are total number of measurements of the time series. The ex- 
pectation value of O is an ordinary number, meanwhile the estimator O is 
a random number fluctuating around the theoretically expected value. For 
numerical simulations the estimation would require to repeat the simulation 
several times, forced by the limit of computing power, the fluctuations of the 
mean value are not probed directly. The usual way is to estimates its variance: 



using the distribution of the individual measurements of O. In case of uncor- 
rected subsequent T MC measurements the previous formula reduces to: 



(B.l) 



i=l 




2 



(B.2) 




(B.3) 
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Oq is the variance of the individual measurements. Explicitly the (1B.3j) ex- 
pression is: 



To formulate the (IB.2|) we assume time-translation invariance over the Markov 
chain and the equilibrium of simulation (values taken after the thermalization) 
for any distribution P(0) of the O. It is usual to refer to oq as the variance 
of the mean, it is the squared width of this distribution and is used as the 
squared error quoted together with the value O. In the case of a Gaussian 
distribution we have a probability of about 68% that all identical simulations 
will produce a mean value in the interval [O — gq, O + erg]. If we use 2cr as the 
error the percentage become about 95% and for 3a interval the probability is 
higher than 99.7%. 

If the samples are not statistically independent flB.3|) becomes an under- 
estimation of the error. The correct expression for erg is [77J: 



*o = \l (O 2 - {Of) (B.5) 

where AT is the Monte Carlo time interval (T = Tmc/AT) and r is called 
autocorrelation time. The autocorrelation time is defined as 1791: 



T 2 



k=i ^ ' 



with 



i=i 



For large T and 2r, equation (1B.5P turns into: 



Vo = \l^(o*-(Oy) (B.8) 



T. 



MC 



The equation (1B.5|) states that due to temporal correlations of the measure- 



ments the real statistical error aQ of the Monte Carlo estimator O is multiplied 
by a factor of \/2t. We can reformulate the correlated statistical error looking 
to the uncorrelated CclSG clSi 




°'6 = \h^ (B.9) 
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but now with a effective statistic parameter 

T eff = < T MC (B.10) 

The effective statistic, for correlated data, is smaller of the number total num- 
ber of measurements and are approximately uncorrelated for 2r iterations. To 
obtain an acceptable error estimation (especially for such quantity like specific 
heat and susceptibility) a simulation has to provide an estimation of autocor- 
relation times too. 



B.l Methods to estimate the error 



B.l.l Binning method 

In many cases the absence of autocorrelations time factor on statistical errors 
bring an great underestimation, beside the calculus of autocorrelation time 
usually involves a big amount of computation time. The statistical error un- 
derestimation can occur even for completely uncorrelated data in time, in fact 
arises the problem of the error estimation for quantities that are not directly 
measured in the simulation but are derived as a general combination of direct 
measures, but in this case can be use the error propagation. An good balance, 
in some situation, can be found using some alternative methods to estimate 
the statistical errors; the binning analysis is much more easy to compute but 
is less accurate of the direct computation. This method consist to split the 
original correlated measurements vector {O^Jl'f into smaller non-overlapping 
T B blocks of length k such that T MC = T B k and compute the block average 
Ob, u of the n-th block 

1 k 

Ob,u = r ^2 0(n-i)k+i (B-ll) 
i=i 

with n = 1, • • • , Tg. The mean value over all block satisfies the identity 
Ob = O. For block length k large enough (k <C r) the blocks become approxi- 
mately uncorrelated, therefore the variance can be computed according to the 
uncorrelated estimator flB.4j) . Resuming, we have for the squared statistical 
error of the mean value: 



a o 



2 _ a B 



2 i T B 



^(0 B , n - Ob) 2 (B.12) 



T B T b (T b - 1) 

v ' n=l 

To obtain an acceptable result we must probe various partition of Tb and 
select the one which maximize the statistical error. The main disadvantage 
of this method is the strongly dependence of the error on the choice of the 
partition T B 
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B.1.2 Jackknife estimation 

This approach [8"4"l I8"5] is a variation of the binning method and can be viewed 
cis db r<3~ sampling. This method too concerns to separate the vector of mea- 
surements into Tb blocks but instead consider small blocks of length k, in 
the Jackknife method we build Tb large Jackknife blocks Oj jn containing all 
measurements but omitting the data contained in the previous n-th binning 
block: 

Oj, n = TM f ~ k ° B > n (B.13) 

J-MC - & 

with n = 1, • • • , Tb- The Jackknife blocks consists of Tmc ~ k data, in this 
way the n-th block contains almost every measurements of the complete vec- 
tor. This allows to use this method even to estimate the errors in non-linear 
combinations of measured variables, the bias is comparable to that of the total 
data. This choice of the Tb Jackknife blocks makes them trivially correlated 
due to the common data coming from the same original data in Tb — 1 dif- 
ferent Jackknife blocks. However, this correlation is a pure re-sampling and 
it is not connected to time correlation. As a consequence, it does not mag- 
nify the error but the Jackknife block variance <jj will be smaller than the 
variance computed using binning method. This reduction can be corrected by 
multiplying <jj with a factor (Tg — l) 2 , obtaining: 

n=l 



B.2 Numerical estimation of autocorrelation 
time 

A third option, the so called Sokal-Madras method [83], is an hybrid approach 
founded on the estimation of the autocorrelation time using (1B.6j) - (lB.7j) . In 
the case that the autocorrelation time satisfies r < 0.5 the samples are not 
correlated thus, to estimate the error, is used the standard uncorrelated error 
given by (1B.4P otherwise error is given by (jB.7[) . 

An estimation of the autocorrelation time is important not only for the 
error estimation, but also for the computation of static quantities themselves. 
In general, it is very hard to have a priori estimation and a direct numerical 
analysis is usually too time consuming. An estimator A(k) for the autocorrela- 
tion function is obtained by replacing in (1B.7P the expectation values by mean 
values. Can be shown that increasing k the variance of A(k) grows quickly, to 
achieve an rough estimation of r and thus a correct error estimate, it is useful 
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to record the partial autocorrelation time estimator 

r\k max ) = -+J2Mk) (B.15) 

K=l 

In the limit of large k max the previous sum tents to r, unfortunately its 
statistical error increases rapidly. An usual procedure, used to find a rea- 
sonable balance between systematic and statistical errors, is to determine a 
self-consistently upper limit k max by cutting off the summation for k max > 
6t' (k max ) . 

In this thesis were compared the errors obtained using the previous three 
methods. Referring to the basic observables like energy density E and the 
order parameters like ifi, Zi^, we have found that the errors were compatible 
among the three methods. For quantities like the specific heat the errors were 
no more compatible, therefore was applied the more accurate Sokal-Madras 
method using the described cut for the estimation of the correlation time. 



B.3 The update algorithm 

The model considered in this works involves five complex independent matrices 
for the 4-dimensional case and 3 complex independent for the 2-dimensional 
one, beside the approximated actions are quite complex and contain terms 
up the fourth power. The usual update algorithm, used to generate a new 
configuration, expects to change all the entries of the independent matrices: 

ipij — > iJij+aoij, Z kij — > Z kij +a kij , with k = 0, • • -4, i,j = 1, • • -N (B.16) 

where N is the matrix size and a k ij are random complex numbers in which the 
real and imaginary parts vary in [—2, 2]. The amplitude of the random inter- 
val was numerically determined to obtain stable results for the observables for 
larger interval. The new proposed configurations are judged by the metropo- 
lis algorithm, which requires to recompute the action for the new proposed 
configuration. Using an optimized algorithm, this operation involves for each 
Metropolis check, about 16iV 3 or 6N 3 (taking in account the cross products 
between the independent matrices present in the actions) as a leading order 
of number of operations respectively for the 4-dimensional and 2-dimensional 
case. In order to reduce the computation time, instead update all entries 
we have changed just one matrix coefficient for each Monte Carlo step, the 
coefficient to update is not chosen randomly but following an order used to 
vary all the coefficients after 5N 2 or 3N 2 steps. In this way to execute the 
Metropolis check we can compute the variation of the action respect the pre- 
vious configuration with a leading order of operations proportional to N. The 
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drawback of this approach are manly two: the complication of the code and 
the behavior of the correlation time in N. In fact two configurations will share 
at least one coefficient until 5N 2 or 3N 2 Monte Carlo steps, this correlations 
introduce an increasing of the correlation time increasing N and presumably 
r behaves like 5N 2 , 3N 2 . However, using this optimization, together with the 
already showed theoretical simplifications and the implementation of parallel 
computing, we were able to execute our simulation for N = 20 with sufficient 
precision and in reasonable time despite our limited computation resource. 
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A' 2 

0.25 1 
0.20 1 
0.15 1 
0.10 1 
0.05 1 



f 
N 

0.30 - 
0.25 ; 

0.20 ; 
0.15 ; 

0.10 1 
0.05 - 



